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Abstract The X-Ray Telescope (XRT) onboard the Hinode satellite is an X-ray 
imager that observes the solar corona with unprecedentedly high angular resolu- 
tion (consistent with its 1" pixel size). XRT has nine X-ray analysis filters with 
different temperature responses. One of the most significant scientific features 
of this telescope is its capability of diagnosing coronal temperatures from less 
than 1 MK to more than 10 MK, which has never been accomplished before. To 
make full use of this capability, accurate calibration of the coronal temperature 
response of XRT is indispensable and is presented in this article. The effect of on- 
orbit contamination is also taken into account in the calibration. On the basis of 
our calibration results, we review the coronal-temperaturc-diagnostic capability 
of XRT. 

Keywords: Corona; Instrumentation and Data Management 



1 National Astronomical Observatory of Japan (NAOJ), 

2- 21-1 Osawa, Mitaka, Tokyo, 181 -8588, Japan 
email: noriyuki.narukage@nao.ac.jp 

2 Institute of Space and Astronautical Science/ 
Japan Aerospace Exploration Agency (ISAS/JAXA), 

3- 1-1 Yoshinodai, Sagamihara, Kanagawa, 229-8510, Japan 

3 Nobeyama Solar Radio Observatory/ 
National Astronomical Observatory of Japan, 
Minamimaki, Minamisaku, Nagano, 384-1305, Japan 

4 Space Engineering Development Co., Ltd., 
1-12-2 Takezono, Tsukuba, Ibaraki, 305-0032, Japan 

5 Smithsonian Astrophysical Observatory, 

60 Garden Street, Cambridge, MA 02138, USA. 

6 NASA/Marshall Space Flight Center, 
AL 35812, USA 



SOLA: Narukage_etal_R3.tex; 15 November 2010; 1:16; p. 1 



N. Narukage et al. 




X-ray telescope (XRT) 



Figure 1. The Hinode satellite has three telescopes, namely Solar Optical Telescope (SOT), 
EUV Imaging Spectrometer (EIS), and X-Ray Telescope (XRT), to observe the Sun in different 
wavelengths. The XRT is mounted at the bottom of this image, as indicated. The size of Hinode 
is approximately 2 m (from XRT to EIS) X 10 m (between both ends of solar paddles) X 4 m 
(between Sun- facing and rear ends). 



1. Introduction 



The outer atmosphere of the Sun, the solar corona, is most clearly discernible 
when seen in soft X-rays. Since the early rocket experiments in the late 1960s, 
it has become widely recognized that soft X-ray imagery of the Sun provides a 
powerful means to investigate physical conditions of hot plasmas (whose tem- 
perature often exceeds 1 MK) that prevail in the corona. Coronal imaging 
with the Soft X-ray Telescope (SXT) onboard the Yohkoh (Solar-A) satellite 
(Og awara et al, 1991||Tsuneta et al., 1991 1 (operation period: 1991-2001) cov- 
ered a full solar activity cycle and have revealed that magnetic reconnection plays 
an essential role in the energy release processes in the dynamic solar corona. 

The SOLAR-B satellite was launched at 21:36 UT on 22 September 2006. 
It was named "Hinode ", which means sunrise in Japanese ( |Kosugi et al.\ 120071 
Figure [lj . The X-Ray Telescope (XRT) onboard Hinode is a successor of the 
Yohkoh /SXT. It also employs grazing- incidence optics, but with improved spatial 
resolution (consistent with 1" CCD pixel size compared with ps 2.5" pixel size 
in the case of SXT) while maintaining similar exposure cadence to that of SXT 
(only order of few milli-seconds for flare and order of few seconds for active 
region) . The XRT has the capability of imaging emission formed at much longer 
wavelengths than SXT. The combination of a backside thinned CCD and thin 
Al-mesh filter allow XRT to image emission significantly longward of 60 A. This 
is a major difference from Yohkoh /SXT in that the XRT can observe not only 
high-temperature plasmas (> 2 MK) seen with SXT, but also low-temperature 
(< 2 MK, reaching even below 1 MK) plasmas which comprise a significant 
amount of the corona. With this enhanced temperature range for observing 
coronal plasmas, coupled with increased spatial resolution, the XRT is able to 
perform detailed imaging observations of a wide variety of coronal plasmas in 
a temperature range covering, continuously, from below 1 MK to well above 
10 MK. One of the most significant scientific features of the XRT is its coronal- 
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temperature-diagnostic capability, namely capability to make temperature maps, 
for such plasmas. 

In order to have XRT perform coronal temperature diagnostics with its full 
capability, we carefully calibrated the effective area of the XRT and its response 
to coronal temperatures using not only ground-based test data but also on-orbit 
data observed in X-rays and visible light. The effect of the on-orbit contami- 
nation, which manifested itself as decreasing intensity of the Sun's corona with 
time as imaged by XRT, was also calibrated as accurately as possible. On the 
basis of our calibration results, we review the coronal-temperature-diagnostic 
capability of XRT with the filter-ratio method. 

In Section [5J the optical elements of XRT are briefly mentioned. In Section [31 
we show the overview of the calibration performed in this article. In Sectional we 
summarize how to identify the contaminant and how to measure its accumulating 
thickness on the focal-plane analysis filters (FPAFs) and CCD as a function 
of time. This then identifies the effective area and temperature response of 
XRT including the contamination at every period. In Section [5J we evaluate 
the coronal-temperature-diagnostic capability of XRT, explain the filter-ratio 
method to derive the coronal temperature, and present suitable filter pairs for 
each coronal temperature. An example of coronal temperature distribution with 
the XRT data is also shown. Finally, we summarize the result of this article in 
Section [51 

Additionally, in Appendix [K\ the calibration of the nine X-ray FPAFs with 
ground-based end-to-end test data is explained. In 2002-2003, the X-ray trans- 
mission measurement of the FPAFs was performed at the X-ray Astronomy 
Calibration and Testing (XACT) facility of the Istituto Nazionale di Astrofisica / 
Osservatorio Astronomico di Palermo "G.S. Vaiana" ( |Collura et ai, 1994 ). How- 
ever, some of the calibrated FPAFs were unfortunately damaged. In this article, 
we characterize the complete flight set of FPAFs, namely both the un-damaged 
and re-manufacturcd FPAFs, with another ground-based test performed in 2005. 
Our calibrated thicknesses arc consistent with the thicknesses available from 
XACT data. 

An overview of the XRT scientific objectives, design, and performance of the 
XRT telescope are summarized bv lGolub et al.\ (|2007|) . while the X-ray camera 
is described bv lKano et~aL\ (|2DD5|> . 
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Figure 2. Optical elements of XRT. The solid line shows the optical path in X-rays. The 
optical elements for X-rays are the pre-filter (PF), primary mirror (Ml), secondary mirror 
(M2), focal-plane analysis filters (FPAF1 and FPAF2) mounted on filter wheels (FW1 and 
FW2, respectively), and CCD. The optical path in visible light (G-band) is indicated by the 
dashed line. The optical elements for visible light are front G-band filter, G-band lens, G-band 
filter mounted on filter wheel 2 (FW2), and CCD. 

2. Optical Elements of XRT 

In order to describe in detail the calibrations performed on the XRT, we first 
briefly describe the optical elements. Figure [2] shows the optical elements and 
optical paths of XRT. The XRT can take not only X-ray images but also visible 
light (G-band) images for co-aligning images between XRT and other instru- 
ments mounted on Hinode. The solid and dashed lines in Figure [5] show the 
optical paths in X-rays and visible light, respectively. 

The XRT optical elements for X-rays consist of pre-filter (PF), two grazing- 
incidence mirrors (Ml and M2), two focal-plane analysis filters (FPAF1 and 
FPAF2), and CCD. 

2.1. Optical Elements for X-rays 

i) Pre-filter: The pre-filter consists of eight fan-shaped annular filter segments 

each made of thin aluminum on a polyimide film (Figure [3]) , which were 
manufactured by Luxel Corporation. The aperture shape was designed to 
obscure the locations of mirror-bonding pads to avoid possible degrada- 
tion in image quality caused by surface deformation of the mirror at the 
bonding pads. In Figure [31 the aperture areas are shown in color with their 
respective opening angles. The average thickness of the entire pre-filter 
assembly is then given as XX ^ x S)/'^2,0, where d and 9 are the thickness 
and opening angle of each filter segment, respectively. The average pre-filter 
consists of 1538 A Al and 2030 A polyimide based on the measurement by 
manufacturer. 

ii) Mirrors: The Ml and M2 mirrors are grazing-incidcncc annular mirror 

manufactured by Goodrich, each made of Zerodur. On the basis of mea- 
surements by the manufacturer, the annular entrance aperture of the XRT 



SOLA: Narukage_etal_R3.tex; 15 November 2010; 1:16; p. 4 



Coronal Temperature Diagnostic Capability of the ffinode/X-Ray Telescope 



Position 3 
Al 1500A 
poly 2044A 
angle 



19.87 



Position 5 
Al 1527A 
poly 2031 A 
angle 34.87 



Position 7 
Al 1527A 
poly 2031 A 
angle 40.94 



Position 1 
Al 1502A 
poly 201 9A 
angle 40.94 



Position 8 
Al 1600A 
poly 1988A 
angle 4.87x2 




Position 2 
Al 1599A 
poly 2035A 
angle 



19.87 



Position 4 
Al 1502A 
poly 201 9A 
angle 34.87 



Position 6 
Al 1601 A 
poly 2050A 
angle 40.94 



Figure 3. Entrance aperture and prc-filtcr of XRT. The pre-filter consists of eight fan-shaped 
annular segments made of Aluminum on a polyimide film. The filter thicknesses measured by 
Luxel and the opening angles are summarized. 



filter wheel 1 filter wheel 2 




Figure 4. Filter wheels and focal-plane analysis filters. XRT has two filter wheels: Filter 
Wheels 1 and 2 (FW1 and FW2). Each filter wheel has six positions to mount the filters. Five 
X-ray analysis niters are mounted on FW1, and FW2 has four X-ray filters and one filter for 
visible-light observation (G-band). The position number and the name of mounted filter are 
described in each filter position. 
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primary mirror is located between radii of r\ = 17.042446 cm and T2 = 
17.074051 cm. Considering the open angle of pre-filter (242.04°, see Fig- 
ure [3]), the geometric aperture area of XRT is calculated to be 7r x (r\ — 
rf) x (242.04/360) = 2.28 cm 2 . The grazing angle of X-rays at each XRT 
mirror is about 0.91° on average. 
Hi) Focal-plane analysis filters: XRT has two filter wheels: FW1 and FW2. 
Each filter wheel has six positions to mount focal-plane analysis filters 
(FPAFs). As shown in Figure 01 five X-ray analysis filters arc mounted 
on FW1, and FW2 has four X-ray filters and one filter for visible-light 
observation (G-band). These nine X-ray filters arc made of several kinds 
of metal and support, with different thicknesses as summarized in Table [TJ 
The filters are designed to observe the corona in a temperature range from 
less than 1 MK to more than 10 MK. Focal-plane analysis filters on FW1 
and FW2 (hereafter FPAF1 and FPAF2, respectively) can be selected inde- 
pendently, even allowing combinations of filters both from FW1 and FW2 
in series. 

iv) CCD: X-rays are focused on a back-illuminated CCD ( |Kano et al., 2008| . 

The effective area [A a s] of XRT is defined by the product of the geometric 
aperture area [^4] and the efficiency of all of the optical elements: 

A e g = A X 7p F X R M 1 X R M 2 X 7FPAF1 X 7FPAF2 x QEqcb, (1) 

where 7pf is the transmission of the prc-filter, Rmi and Ryi2 reflectivities at the 
primary and secondary mirrors, 7fpafi and 7fpaf2 the transmissions of FPAF1 
and FPAF2, and QEqcu the Quantum Efficiency of the CCD. Because each 
efficiency is a function of wavelength, A e g is also a function of wavelength. In 
this article, we calibrate Tpf, Tfpafi, and 7fpaf2, and then derive the effective 
area. 

2.2. Optical Elements for G-band 

When XRT takes visible-light images in G-band, the visible-light shutter is 
opened, and FW1 and FW2 are set to the open and G-band filter position, 
respectively This G-band filter is a bandpass filter with a central wavelength 
of 4305.6 A and a bandwidth of 172.8 A (FWHM). The visible light is focused 
onto the same CCD as X-rays (Figure [2]). X-rays cannot reach the CCD through 
the G-band (dotted line) or the X-ray (solid line) paths, because those X-rays 
are blocked by G-band filters employed at the front of the XRT and on FW2, 
respectively. 
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Table 1. Focal- plane analysis filters and pre-filter. 



FW-posM 


filter name 


pre-delivery measurements at Luxel'- 6 ' 




calibrated values in this article 




metal 


support 


pure metal 


oxidized 


at fabrication 1 -^ 


support 


1-0 


open 














1-1 


Al-poly 


Al< d ) 1283A (±50A) 


poly( e ) 2656A (±100A) 


Al 1412A 


AI2O3 75A 


1470A 


poly 2656A 


1-2 


C-poly 


CC/) 6038A (±50A) 


poly 3478A (±100A) 


C 5190A 




5190A 


poly 3478A 


1-3 


thin-Bc 


BeCs) 9/1 (+5/i) 




Be 10.46/1 


BeO 150A 


10.47/1 




1-4 


mcd-Bc 


Be 30/i (+5/i/-2/i) 




Be 26.89/t 


BeO 150A 


26.90/1 




1-5 


med-Al 


Al 12.5/1 (±5%) 




Al 12.25/1 


AI2O3 150A 


12.26/1 




2-0 


open 














2-1 


Al-mesh 


Al 1605A (±50A) 


meshCO 


Al 1583A 


AI2O3 150A 


1700 A 


77% trans, mesh 


2-2 


Ti-poly 


TiC*) 2345A (±50A) 


poly 2522A (±100A) 


Ti 2338A 


Ti0 2 75A 


2380A 


poly 2522A 


2-3 


G-band 


Glass 












2-1 


thick-Al 


Al 25m (±10%) 




Al 26.09/1 


AI2O3 150A 


26.1/1 




2-5 


thick-Be 


Be 300/i (±30%) 




Be 252.79/t 


BeO 150A 


252.8/1 






pre-filter 


Al 1538A (±50A) 


poly( e ) 2030A (±100A) 


Al 1492A 


AI2O3 75A 


1550A 


poly 2030A 



( a ) "FW" and "pos" mean the filter wheel number and position on the filter wheel, respectively. 

CO These values are described in the certification sheet by Luxel. Note that the thicknesses measured by Luxcl are only for reference. 

O The expected metal thickness when it was fabricated. At fabrication, the metals had not oxidized at all. This value is derived from the calibrated 
thicknesses of pure and oxidized metal using Equation J56D in Appendix U This is shown for comparison with the value measured by Luxel. 

(<0— CsJiW "Al" means Aluminum, "poly" polyimide (C2 2 HioN205), "C" Carbon, "Be" Beryllium, and "Ti" Titanium. 

CO The mesh for Al-mcsh filter is made of stainless steel. 
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3. Overview of the Calibration 

In Section^] and Appendices, we will carefully calibrate the spectral response of 
XRT including the effect of the on-orbit contamination as accurately as possi- 
ble. The on-orbit contamination makes our calibrations complicated. Therefore, 
before proceeding in detail, we first give an outline of calibration activities in 
Figure [5] 

In this figure, time passes from left to right. After the launch of Hinode, from 
the viewpoint of the XRT thermal environment, there are three distinct intervals: 
Phase 1 (22 September 2006 - 19 October 2006), Phase 2 (19 October 2006 - 18 
June 2007), and Phase 3 (18 June 2007 - present). Phase 1 is the period from 
the launch to first light, when the CCD bakeout heater was kept on. Phases 2 
and 3 are the periods of normal operation of XRT without and with enabling 
(turning on) of the operational heater, respectively. The details of each phase 
are described in Appcndix lB.il 

At the bottom of Figure the frequency of observations in X-rays and the 
G-band are shown by fines. The solid lines indicate that many data sets were 
taken, dotted line a few data sets, and no line no signifies data sets. From 
the viewpoint of this XRT observation, Phase 2 is divided into two periods, 
namely the early period (19 October 2006 - 16 February 2007) and the late 
period (16 February 2007 - 18 June 2007). In the early period, there are no 
simultaneous observations in X-rays and G-band, while in the late period, there 
are simultaneous observations. 

The XRT calibrations were performed in order from top to bottom of the 
figure. The white boxes indicate the data sets used for our calibration. The black 
arrows connect the data sets or models to the calibrated results. The calibrated 
XRT instrumental parameters are shown by green boxes. The bars in light-green 
are extended from the green boxes to the times at which the calibrated results 
in the green boxes are valid. Note that the results shown by green boxes and 
light-green bars are used for later calibration without indicating by black arrows. 
The blue and red boxes are the calibration results of contaminant accumulated 
on the CCD and FPAFs, respectively. 

In Phase 3, since both X-ray and G-band data were taken sufficiently often 
(see bottom part of Figure [5]) , the calibrations of both the XRT instrument and 
contamination are possible. However, in Phase 2 the data in X-rays and/or G- 
band were lacking for the calibration. In order to calibrate XRT for such episodes, 
we establish three models shown in yellow boxes. For the late period of Phase 2 
where G-band data are sparse, the model of the "XSC (X-ray-suggested CCD 
contamination) relation" based on the result of Phase 3 is used to calibrate the 
contaminant on both CCD and FPAFs as shown by a bar in light-yellow. Mean- 
while, for the early period of Phase 2 where the data are sparser than in the late 
period, two models — "ODC (operational heater driven contamination)" and 
"USF (FPAF-usage-suggestcd FPAF contamination)" relations — established 
in the late period are applied. In Phase 1, we expect that the contaminant did 
not accumulate on both CCD and FPAFs, because the CCD bake heater was 
kept on and the temperature around the FPAFs was warm. 

As shown in Figure \5\ we calibrate the spectral response of XRT throughout 
its whole operation period. Table Q] summarizes the calibrated instrumental 
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Figure 5. Outline of calibration activity for XRT. 
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Figure 5. Continued. 



parameters of XRT with ground-based end-to-end test data (shown by green 
boxes in Figure [5]) . The details of this ground-based calibration are described 
in Appendix [X] Time- varying thicknesses of on-orbit contaminant accumulated 
on the CCD (blue boxes) and FPAFs (red boxes) are given by Figures [H and 
[3 respectively. The contaminant material is identified to have the chemical 
composition of a long-chain organic compound without silicon, refractive index 
of ss 1.5, and density of ~ 1 g cm -3 . In Section @] and Appcndix[B] the on-orbit 
calibrations arc described. 
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Figure 6. Calibrated thickness of contaminant accumulated on the CCD, with bakcouts 
indicated by gray vertical lines. 
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Figure 7. Calibrated thicknesses of contaminant accumulated on the FPAFs. The final 
thickness of contaminant on each FPAF is shown. 
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4. On-orbit Calibration 

XRT took its first image on 23 October 2006. After this first light, XRT has 
been taking several thousand images a day. However, several months into the 
mission, we found that the X-ray intensity seen with XRT started to decrease 
continuously, especially, in the thinner filters. 

The top panel of Figure [5] shows the intensity [DN sec -1 pixel -1 ] of the quiet 
Sun in the four X-ray filters of XRT. In an intensity distribution (histogram) 
from a full-Sun corona, quiet-Sun areas manifest themselves as a concentration 
in the distribution whose profile is well expressed by a Gaussian. We took 
the center position of the Gaussian profile as representing the area-averaged 
quiet-Sun intensity. The orange, red, green, and blue pluses indicate the quiet- 
Sun intensity observed with the Al-mesh, Al-poly, C-poly, and Ti-poly filters, 
respectively. The gray area shows episodes of CCD bakeout (sec Table [7] in 
Appendix [Q for details). The intensity decrease is seen in all filters, with that 
in the Al-mesh filter being most significant. While intensity fluctuations in the 
short-term (several weeks) is caused by solar activity such as the appearance of 
active regions, the systematic intensity decrease can not be explained by solar 
activity. For example, the intensity of the quiet Sun with the Al-mesh, Al-poly, 
C-poly, and Ti-poly filters decreases to 27%, 59%, 80%, and 73%, respectively, 
during 142 days from 8 September 2007 to 28 January 2008. We conclude that the 
decrease is caused by accumulation of materials, namely contaminants, obscuring 
the optical path of the XRT. 

The X-ray intensity ratio is a function of the XRT filter response and coronal 
temperature (see Equation ©), and is not affected by the variation of the coronal 
density nor the distance between the Sun and Hinode unlike the intensity. Since 
such intensity ratio recovered to almost the same level by each CCD bakeout 
as seen in Phase 3 in the bottom panel of Figure [51 it is most likely that at 
least some fraction of the contaminants is accumulating on the CCD and each 
bakeout reduces the contaminant to the same thickness. The analysis of G-band 
intensity in Appendix [B] gives the result that this "same thickness" is actually 
zero. Hence, we can say that each bakeout completely removes the contaminant 
from the CCD (except for the "spots" discussed in Appendix IB"]) . Details of CCD 
bakeout are given in Appendix [Cl 

We also find the intensity of the quiet Sun taken with Al-poly filter (red 
pluses) is weaker than with the C-poly filter (green pluses) as clearly seen in 
the September 2007 to February 2008 time frame in the top panel of Figure HJ 
Strangely enough, however, the ground-based measurement at XRCF indicates 
that Al-poly filter should have a larger signal over the C-poly filter throughout 
all coronal temperatures (see the top panel of Figure fT0|) . This discrepancy led 
us to the notion that, in addition to the CCD, the FPAFs also suffered from the 
accumulation of contaminants with different accumulated thickness for different 
filters. 

Here, we briefly summarize how we characterized the material of the contam- 
inant and the accumulation thicknesses on the CCD and FPAFs as a function of 
time. First, we measure the thickness of contaminant accumulated on the CCD 
using G-band data and unique method (see Appendix IB . 2 1 in detail) as shown in 
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Figure 8. Intensity of the quiet Sun observed with XRT. top panel: Original intensity of 
the quiet Sun observed with four kinds of X-ray filters, bottom panel: Intensity ratio of the 
quiet Sun normalized with the intensity observed with Ti-poly filter. The gray area shows the 
duration of CCD bakcout. 

Note - The information about the phase, which is defined in Appcndix lB.il and status of bake 
out heater and operational heater (gray bars mean that these heaters are on) are shown in the 
top part of these plots. Hinode has been placed in a Sun-synchronous polar orbit that is above 
the day-night line on Earth, where Hinode can continuously observe the Sun for 24 hours a day 
for about nine months per year. The remaining three months are called the "eclipse season", 
where the Sun is eclipsed by Earth for a maximum of 20 minutes in each 98-minute orbit. This 
eclipse season is indicated by light-gray bars. 
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Table 2. Behavior of contamination 



period'"' 


bake 


operational 


eclipse 


contamination' 6 ' 




heater 


heater 


season 


on CCD 


on FPAFs 


Phase 1 


ON 


OFF 


no 


none 


none 


normal Phase 2 


OFF 


OFF 


N/A 


/ 




CCP in Phase 2 


OFF 


ON-s-OFF 


no 


//V 


/// 


WCP in Phase 2 


OFF 


ON-s-OFF 


close to 


/W 




Phase 3 


OFF 


ON 


N/A 


/W 




CCD bakeout 


ON 


N/A 


N/A 


completely removed 





( a ) Phases 1-3 are defined in Appendix IB. II CCP and WCP represent "cool contamination 
periods" and "warm contamination period" defined in Appendix IB, 3. 3l respectively. 

(f>) ii y*. yi yv> an( j iiym indicate the rapid and slow accumulation of contaminant, respectively. 
" — >" means that the thickness of contaminant is stable. 



Figure [S] Next, on the basis of the temporal evolution of the G-band and X-ray 
intensities, we identify material of contaminant as the chemical composition of 
a long-chain organic compound without silicon, refractive index of ~ 1.5, and 
density of « 1 g cm -3 . Thickness of the prc-filter, oxidization thickness of the 
pre-filter and FPAFs, and thickness of contaminants accumulated on each FPAF 
are also calibrated (see Appendix IB. 31 in detail) as summarized in Table [1] and 
Figure [7] Table [5] shows the summary of the on-orbit contamination behavior. 

Since we now have accurate knowledge of the response parameters (filter 
thicknesses, contaminant thicknesses, etc.) at every epoch after launch of Hinode, 
we can accurately derive the XRT effective area and the temperature response. 
(The temperature response is explained in detail in Appendix ID.2I ) The cali- 
brated effective area and temperature response are shown in Figures [9] and 1101 
respectively. Let us next look into the temperature-diagnostics capability of XRT 
with the calibrated temperature response. 

5. Temperature Diagnostics with Filter-Ratio Method 

The coronal-tempcrature-diagnostic capability is the most significant scientific 
attribute of XRT. In Scction[57TJ we explain the coronal temperature diagnostics 
with the filter-ratio method and the estimate of statistical errors in the derived 
temperature due to photon noise for the XRT. The filter-ratio method derives a 
mean temperature weighted by the filter responses and emission measure from 
images observed with t wo different filters (|Vaiana, Krieger, and Timothy, 1973[ ) . 
IHara et aZ.I (|1992|) and lKano and Tsunetal (|1995[) discussed coronal temperature 
diagnostics and its statistical errors for the Soft X-ray Telescope (SXT) onboard 
Yohkoh (SOLAR- A: \Tsxmeta e£ "oil [TMTj) . respectively In Section ICT we sum- 
marize the suitable filter pairs of XRT for coronal temperature diagnostics with 
the filter-ratio method. In Section 15.31 we discuss the meaning of filter-ratio 
temperature derived with XRT. 
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Figure 9. Effective area of XRT on orbit considering the contaminant accumulated on the 
FPAFs and CCD. The solid and dotted lines show the effective area just after CCD bakeout and 
one month after the bakeout in Phase 3, respectively. We assume that 800 A of contaminant 
accumulated on the CCD one month after the bakeout. 



5.1. Filter- Ratio Method 

Assuming that the corona is isothermal, the data number (DN) observed with 
XRT is 

FIT) 

DN = ——t x VEM x t (2) 

for an isothermal corona at a temperature of T, where F(T), VEM, t, and 5 are 
the XRT temperature responses, volume emission measure, exposure time, and 
solar area detected in 1 pixel of CCD (ss 560000 km 2 derived from plate scale 



of XRT ps 1.03 arcsec pixel -1 ( Shimizu et al, 2007 ) and the relation of 1 arcsec 
ps 726 km arcsec -1 on the solar surface), respectively. Then, the ratio i?y (T) 
of normalized DNs per unit time taken with two different filters is a function of 
only the coronal temperature T: 



_ DNj/U = Fj(T)/S x VEM F^T) 
>- DNj jtj Fj (T)/Sx VEM F (T) ' 



The subscripts i and j specify the two different filters. 

The coronal temperature can be estimated with the observed ratio of normal- 
ized DNs as 
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Figure 10. top panel: Response of XRT to the coronal temperature at the launch of Hin- 
ode, i.e. when XRT had not been contaminated, bottom panel: Response of XRT to coronal 
temperature on orbit considering the contaminant accumulated on the FPAFs and CCD. The 
solid and dotted lines in the bottom panel show the temperature response just after CCD 
bakeout and one month after the bakeout in Phase 3, respectively. We assume that 800 A of 
contaminant accumulated on the CCD one month after the bakeout. 
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temperature volume emission measure 




1.0 1.5 2.0 2.5 [MK] 41 42 43 44 [log cm 3 ] 



Figure 11. Example of the XRT temperature diagnostics with the filter-ratio method applied 
to full-disk images taken with Al-mesh and Al-poly filters on 12 February 2007. Top-left panel: 
the original X-ray data taken with Al-poly filter. Top-right panel: the ratio of temperature 
response between Al-poly and Al-mesh filters. Bottom-left panel: the full-Sun temperature. 
Bottom-right panel: the volume emission measure. 



where is the inverse of Rij. Once the coronal temperature is obtained, we 
can derive the volume emission measure as 

VBMmJwW-5g&c.2ML. (5) 
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where n e , Tin, and dl are the electron number density [cm ], hydrogen num- 
ber density [cm -3 ], and unit length along the line-of-sight [cm], respectively. 
Figure [TT] is an example of coronal temperature diagnostics with the filter- 
ratio method applied to XRT data taken with Al-mesh and Al-poly filters on 
12 February 2007. The bottom-left and bottom-right panels are the full-Sun 
temperature and emission measure maps, respectively. 

Next, we discuss the estimate of statistical errors in the derived temperature 
and emission measure due to Poisson noise of incident photons into XRT. Since 
Poisson noise of incident photons dominants over other noise sources (e.g., the 
camera system noise and dark current noise reported bv lKano et al.\ (|2008p ). we 
consider only Poisson noise for the error estimate. Note that lKano and Tsunetal 
(1995) estimated those errors for Yohkoh/SXT. The differences between our 
estimate and that of IKano and Tsunetal (|1995p are the adopted solar spec- 
trum database and its units. We adopt solar photon-number spectra in units 
of [photon cm -2 sec -1 sr - 1 A -1 ] fr om the CHIANTI atomic database version 
6.0.1 (jDere et alX 119971 120091 which is the latest version, when we analyzed) 



with ionization equilibrium from chianti.ioneq (Dere, 20071 and abundance from 



sun.coronaLext.abund (Fcl dman et al, 19 92; Landi, Fcld man, and Dere, 2002[ ), while| 
IKano and Tsunetal ()1995[) used solar emissivity [erg cm -2 sec -1 sr -1 A -1 ] from 

Mewe's spectral data flMewe, Gronenschild, and van den Oord, 1985; M cwe, Lemen, and van den Oord, 1986 

First, we derive two conversion factors as preparation for the error estimate. 
When XRT observes an isothermal corona at a temperature of T and a volume 
emission measure of VEM with p pixels of the CCD and t [sec] exposure, the 
incident photon number in a range of wavelength from A to A + dX is 

dN = (^P iso0 (A, T) x s x ^£J^1 x dxj x | x VEM, (6) 

where Pi SO (A,T) [cm™ 2 sec -1 sr -1 A -1 ] is the photon-number spectrum emit- 
ted from an isothermal corona at a temperature of T and an unit of column 
emission measure (CEM) of 1 cm™ 5 , s is the area of 1 CCD pixel (13.5 2 ^m 2 , 
see IKano et al\ l|2008jl ). A c g (A) is the effective area of XRT, and / is the focal 
length of XRT (2708 mm, see IGolub e~taT\ <|2"00"7|l ). The total photon number 
observed with XRT is, then, 



N = J dN. (7) 
The data number (DN) observed with XRT is derived as 

x - x VEM (8) 



(from Equations @ and in Appendix ID. 2[) . where h, c, e, and G are Planck's 
constant, speed of light, elementary electric charge, and the system gain of the 
CCD camera, which is 57.5 [electron DN -1 ] in the XRT case ( |Kano et al, 2008[ ), 



SOLA: Narukage_etal_R3.tex; 15 November 2010; 1:16; p. 18 



Coronal Temperature Diagnostic Capability of the ffinode/X-Ray Telescope 



conversion factor K (1) on orbit 
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temperature [K] 
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Figure 12. top panel: Conversion factor [K^ (T)] of each FPAF from observed data number 
(DN) to incident photon number [N], namely TV = DN/K^ (T). bottom panel: Conversion 
factor [K^ (T)] of each FPAF from observed DN to the square of DN error [ox>jv], namely 
a DN = K(2) ( T ) X DN - These fi 

gures show the conversion factors on orbit considering the 
contaminant accumulated on the FPAFs and CCD, where the solid and dotted lines indicate 
the conversion factors just after CCD bakcout and one month after the bakeout in Phase 3, 
respectively. We assume that 800 A of contaminant accumulated on the CCD one month after 
the bakcout. 
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respectively. When we define K^> (T) as 

/ (p iso0 (A,T) x A cS (A) x j x ^^c) 



( T ) 



dA 



/(p iso0 (A,r)x4 ff (A)) dA 



(9) 



we can rewrite Equation (0 as 



N 



DN 



KM (T) ' 



(10) 



(j 1 ) i s a conversion factor from observed DN to incident photon number. 
The top panel of Figure ^indicates (T) for each filter. 

Since diV follows a Poisson distribution, the standard deviation da of dN is 
expressed as 



dcr 



Since Equation ([5J is rewritten as 



DN 



'dN. 



1 



A ex 3.65 x G 
we can derive the error <jdn of DiV as 



dN, 



(11) 



(12) 



&DN 



he 1 

T X e x 3.65 x G 



) (d*) 2 = y/KW (T) DN, (13) 



where 



^(2) ( T ) 



P iso0 (A,T)x4 ff (A)x(fx^ 



dA 



/(P iso0 (A,T)x4 J ff(A)x^x i ^< ? 



dA 



(14) 



if ( 2 ) (7 1 ) i s a conversion factor from observed DN to the square of the DN error. 
The bottom panel of Figure [T2l shows (T) for each filter. 

We note that these two conversion factors K^ 1 ' and vary in a part 

of detectable temperature range of XRT from ss 1 MK to « 5 MK (sec Fig- 
ure IT2"]) . while those factors for Yohkoh /SXT was almost constant at three 
throughout the whole detectable temperature range of w 2 MK to > 10 MK 
( |Kano and Tsuneta, 1995[ ). 

Following the Appendix of IKano and Tsunetal (| 1995[) , we can estimate the 
error of derived temperature and volume emission measure from these conversion 
factors as 



~T 



dinner) 
dlnT 



(2) 



\ DN 



DN, 



(15) 
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Table 3. Typical values of coronal structures 



coronal 


DEM 


electron 


exposure 


binning 


structure 


reference 


number density 


time 


size 


flare 


IDere and Cookl 119791 


10 110 cm" 3 


1 sec 


1 x 1 


active region 


IWarren et al.\ (120101 


10 9 - 4 cm" 3 


10 sec 


2x2 


quiet Sun 


IBrooks and Warrenl 112006! 


10 8 - 3 cm" 3 


30 sec 


4x4 


coronal hole 


IVernazza and Hccvcsl (|1978|l 


10 8 cm" 3 


64 sec 


8x8 



and 



Qvem 
VEM 



dlnfljj (T) 
dlnT 



\ 



" dln^(T) - 
dlnT 



k[ 2) (T) 

DN, 



'dlnFj(T)' 
dlnT 



Kf (T) 
DNi ' 



(16) 

These relations show how the errors in observed data number [DN] propagate to 
the errors in the derived temperature [T] and volume emission measure [Vi?Af] . 

We note that the temperature in the right-hand side of Equations (fT5|) and 
(fT6l) should be the actual temperature of the coronal plasma. However, we get the 
estimated coronal temperature with the filter-ratio method. Hence, we used the 
derived temperature to estimate the errors instead of the actual temperature. In 
the case of poor signal-to- noise ratio (S/N) data, the difference between actual 
and derived temperatures may be larger than in the case of good S /N data, and 
the above treatment of T in error estimate will cause the additional error. We 
confirmed that 10% S/N data makes this additional error sufficiently small. 

On the basis of the error estimate discussed above, the full-Sun temperature 
map in Figure [TT] is created. In this map, the poor S/N pixels, namely pixels in 
the coronal hole and quiet Sun, are summed (8x8 pixels binning maximum) 
to improve the S/N. Using this good S/N data, the coronal temperature with 
the filter-ratio method can be derived with less than 20% errors across the Sun. 
This diagnostic capability is the most significant scientific feature of XRT. 



5.2. Suitable Filter Pairs for Coronal Temperature Diagnostics with the 
Filter-Ratio Method 



On the basis of the discussion in the previous section, we present suitable filter 
pairs to derive temperatures with filter-ratio method for various coronal struc- 
tures: flares, active regions, quiet Sun, and coronal holes. The typical values for 
each coronal structure are summarized in Table [3] The electron number density 
is calculated as the integral over temperature of DEM with the assumption 
that the line-of-sight depth is 10 5 km. The exposure time and binning size is 
determined from the time scale and spatial scale of each structure, respectively. 
Figures Q2] and HH show the "suitable filter pairs" and their "useful temperature 
ranges" for some representative cases of contamination accumulated on the CCD 
and the FPAFs. The "useful temperature range" means the range where a filter 
pair can estimate the temperature (within 20% accuracy) with the filter-ratio 
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suitable filter pairs just after CCD bakeout 



med-AI/open & 


open/thick-Be 
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Figure 13. Suitable filter pairs to derive the temperature with the filter-ratio method and 
their useful temperature ranges just after the CCD bakeout in Phase 3, namely there is no 
contamination on the CCD but the contaminant accumulated on the FPAFs. The red, orange, 
green, and blue bars show the useful temperature ranges for flare, active region, quiet Sun, and 
coronal hole, respectively. The "v" marks indicate the filter-ratio temperatures derived from 
the DEMs of each region. 

method for a given coronal structure with the typical exposure time and binning 
size (see Table [3]). The "suitable filter pair" means that the pair has a wide 
"useful temperature range" as shown in Figures [Jj|] and For example, the 
Al-mesh and Ti-poly filter pair can have a wide temperature coverage for a wide 
variety of coronal structures, coronal holes, quiet Sun, and active regions. Thus, 
this pair is suitable for full-disk synoptic observations. 

Let us now examine how to derive the suitable filter pairs and their useful 
temperature ranges. Figures IT51 and [TBI are examples of the filter-pair study for 
quiet Sun and active region, respectively A summary of the study is described 
in the top-left of these figures. Here [target] indicates the coronal structure to 
be analyzed. The typical differential emission measure (DEM) of such a target 
is plotted in the top-right side of these figures (the references of typical DEMs 
are summarized in Table [3]) . The row with [binning size] is the recommended 
binning size to collect enough photons. For a bright flare, [binning size] is 1 x 1 
binning, that is no binning. Meanwhile, for the dark coronal hole, [binning size] 
is 8 x 8 binning. The binning size used in our calculation is shown in Table [3] The 
rows with [filter 1] and [filter2] summarize the evaluated filter pair. Labels "fwl" 
and "fw2" indicate used FPAFs on filter wheels 1 and 2, respectively. The label 
"contam." gives the the accumulated thickness of contaminant on each FPAF. 
The label "full-well exp. time" expresses the exposure time necessary to obtain 
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suitable filter pairs one month after CCD bakeout 

med-AI/open & open/thick-Be 
med-Be/open & open/thick-Be 
thin-Be/open & med-Be/open 
open/Ti-poly & med-AI/open 
C-poly/open & med-AI/open 
Al-poly/open & med-AI/open 
open/AI-mesh & med-AI/open 
open/Ti-poly & med-Be/open 
C-poly/open & med-Be/open 
Al-poly/open & med-Be/open 
open/AI-mesh & med-Be/open 
open/Ti-poly & thin-Be/open 
C-poly/open & thin-Be/open 
Al-poly/open & thin-Be/open 
open/AI-mesh & thin-Be/open 
open/AI-mesh & open/Ti-poly 
open/AI-mesh & C-poly/open 
open/AI-mesh & Al-poly/open 

10 5 10 s 10 7 10 B 

temperature [K] 

Figure 14. Suitable filter pairs to derive the temperature with the filter-ratio method and 
their useful temperature ranges one month after the CCD bakeout in Phase 3, namely there 
is contamination not only on the FPAFs but also on the CCD. We assume that 800 A of 
contaminant accumulated on the CCD. The others are the same as in Figure [T3l 

the full- well intensity of the XRT CCD, about 200000 electrons w 3500 DN, when 
XRT observes the coronal structure shown in [target] with [filter 1] or [filtcr2] . The 
full-well intensity means the maximum intensity to keep the linearity between 
the incident photon energy and output intensity as shown by Equation (|33l) . 
In the case of the intensity exceeding the full-well intensity, the actual output 
intensity becomes smaller than the intensity expected from Equation (|33[) . The 
row with [contain, on CCD] indicates the accumulated thickness of contaminant 
on the CCD. 

Next, we explain six conditions that a filter pair shown in [filterl] and [filter2] 
is suitable for the temperature diagnostics with the filter-ratio method. The six 
conditions, (A) -(F), are given below. The four panels at the bottom part of 
Figures [13] and [T|)] show the temperature ranges (white areas) which satisfy the 
following Conditions (A) -(D). 

(A) Expected DN rate at each temperature: 

The middle- left panel shows the expected DN rates (whose unit is DN sec -1 
pixel -1 K _1 ) at each temperature when XRT with [filterl] and [filtcr2] 
observe the [target] which has the DEM shown in the top-right side of the 
figures. Since the dynamic range of the XRT CCD is about 1000, we assume 
that XRT can detect temperature components where DN rate exceeding 
1/1000 of the peak DN rate in each expected DN rate plot. Condition (A) 
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[target] 
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contam. 
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Figure 15. Evaluation of the Al-mesh and Ti-poly filter pair for coronal temperature diag- 
nostics with the filter-ratio method. This is the case of the quiet Sun just after CCD bakeout 
in Phase 3. This pair is suitable for temperature diagnostics in the quiet Sun. 



is that the expected DN rate exceeds the 1/1000 of the peak DN rate (shown 
by dotted lines) in both cases of [filterl] and [filtcr2]. 
(B) Low photon noise: 

The bottom-left panel presents the photon noise, when XRT observes an 
isothermal plasma at each temperature given by the horizontal axis of the 
panel. The electron number density at each temperature is set to be a 
constant value as the typical value described in Table [3] and vertical axis. 
Hence, this plot shows how much photon noise is expected in the observed 
data when we obtain the filter-ratio temperature from such data. The pho- 
ton noise can be calculated with Equation (|13|) as a function of temperature 
and DN. The DN is derived from temperature, VEM, and exposure time 
as shown in Equation ©. For the calculation of VEM, we use the above 
electron number density and 10 5 km as the line-of-sight depth. The effect 
of binning is also considered. The exposure time adopted for this evaluation 
is shorter of the typical exposure time (see Table [3]) or the exposure time 
to obtain the full- well intensity (see the label "full- well exp. time"). The 
adopted exposure time for each filter are described in each filter color in the 
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[target] 
[binning size] 

[filter!] fw1/fw2 
contam. 
full-well exp. time 

[filter2] fw1 / fw2 
contam. 
full-well exp. time 
[contam. on CCD] 
[estimated filter-ratio temperature] 
[useful temperature range] 



ACTIVE REGION 

2x2 binning 
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500 [A] / [A] 
720.75 [ms] 
0[A] 

4.1 x 10 6 
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temperature [K] 
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Figure 16. Evaluation of the Al-poly and C-poly filter pair for coronal temperature diagnos- 
tics with the filter-ratio method. This is the case of the active region just after CCD bakcout 
in Phase 3. This pair is not suitable for temperature diagnostics in the active region. 



bottom-left panel. The horizontal dotted line corresponds to 10% photon 
noise (S/N = 10). We consider that the good S/N data (S/N > 10) should 
be used for the temperature diagnostics, because poor S/N data makes the 
error large as mentioned above. Condition (B) is that the photon noise 
should be lower than 10%. We note that Condition (B) is not so strict, 
because we can sum many images to adequate enough photons. When 
we collect enough photons with a single image, Condition (B) should be 
satisfied. 

(C) Profile of temperature response ratio: 

The middle-right panel indicates the temperature response ratio of [filtcrl] 
and [filter2] shown by Equation ©. Using this ratio, we can estimate 
the coronal temperature from the observed intensity ratio. The horizontal 
dashed lines in Figures [T"5l and \W\ show the expected intensity ratio, when 
XRT observes [target] with [filtcrl] and [filtcr2]. In the case of Figure [T5l 
the temperature is uniquely derived as shown by the vertical dashed line. 
The derived temperature is described by [estimated filter-ratio temperature] 
in the top-left panel of the figure. However, in the case of Figure [TU the 
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temperature cannot be uniquely determined. In this case, [estimated filter- 
ratio temperature] shows the closest temperature to the peak temperature 
of DEM in the top-right panel as a reference. The temperature range where 
the intensity ratio gives the alternative temperatures is indicated by the 
magenta line and gray area. Condition (C) is that the temperature is derived 
uniquely from the filter-ratio method. 

(D) Error in estimated temperature: 

The bottom-right panel shows the fractional error in estimated tempera- 
ture caused by 10% photon noise. The error in estimated temperature is 
calculated with Equation (|T5|) as a function of temperature and DN. The 
DN is derived with the method described in Condition (B). We take an 
estimate error of 20%, which is twice the value of the photon noise, as the 
threshold for a reasonable temperature estimate. This is Condition (D). The 
horizontal dotted line is located at the 20% error in estimated temperature. 

(E) No saturation: 

The shortest exposure time for XRT is 1 ms. Hence, the "full-well exp. time" 
should be longer than 1 ms. For bright coronal structures, namely active 
regions and flares, the suitable filter pairs arc limited by this condition. 

(F) Filter-ratio temperature: 

The useful temperature range with the filter-ratio method is derived from 
satisfying the Conditions (A) - (E) as shown by [useful temperature range] 
and in white area in the top-right DEM panel. We call the filter pair 
suitable if the [estimated filter-ratio temperature] falls within this [useful 
temperature range]. This is Condition (F). 

The suitable filter pairs for coronal temperature diagnostics with filter-ratio 
method should satisfy all of the Conditions (A) - (F) . These suitable filter pairs 
depend on the target, filter pair, and accumulation of contaminant. Examples are 
summarized in Figures [T3] and [H] for some representative cases of contamination 
accumulated on the CCD and the FPAFs with their [useful temperature range] 
and [estimated filter-ratio temperature] shown by bars and "v" marks, respec- 
tively. These suitable filter pairs and their useful temperature range are useful 
for understanding the temperature diagnostic capability of XRT and making 
the XRT observation plan for the temperature analysis. In the analysis with 
actual data taken with XRT, it is also possible to estimate the temperature 
with filter pairs other than those we described here as suitable filter pairs, since 
the suitable filter pairs are derived with the typical values of coronal structures 
(DEMs, density and binning size) summarized in Tabic [31 

We note that since the actual solar corona may have a multi-temperature 
structures, the photon noise for Condition (B) should be estimated with a multi- 
temperature (DEM). However, we consider that the treatment in Condition (B) 
is valid, since we confirmed that the photon noise derived with the filter-ratio 
temperature well represents one derived with DEM, at least, at the estimated 
filter-ratio temperature as will be discussed in Appendix [jj In Condition (D), 
we evaluated the error in estimated filter-ratio temperature. Although the error 
is derived under the assumption that the corona is the single filter-ratio temper- 
ature, the derived error is useful as the accuracy-index of estimated filter-ratio 
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Figure 17. Correlation between DEM-weighted mean temperature [T moan ] and XRT fil- 
ter-ratio temperature [T rat j ] for four DEMs shown in Table [3] This plot is the case when 1 
month passed after the CCD bakcout, i.e. 800 A of contaminant on the CCD. The dashed and 
dotted lines indicate the positions where T rat i is equal to X mea n, and where the difference 
between T mean and T rat i D is 50%, respectively. The number of analyzed filter pairs is shown 
in parenthesis. 



temperature which depends on the photon noise and the shape of temperature 
responses in each filter. 

5.3. Meaning of Filter-Ratio Temperature Derived with XRT 

The actual solar corona may have a multi-temperature structures, and the tem- 
perature derived with the filter-ratio method is a mean temperature weighted 
by the filter responses. Unlike filter-ratio temperatures from narrow-band instru- 
ments (see discussion in |Martens, Cirtain, and Sc hmclz (2002)), filter-ratio tem- 
peratures from broad-band instruments (e.g., Hinode/XRT and Yohkoh/SXT) 
are expected to provide a temperature which is close to the mean temperature 
[Tmean = ^ j ^DE^di? ] weighted by differential emission measure (DEM) of 
the coronal observation target as demonstrated by |Acton, Weston, and Bruner 
for Yohkoh/SXT. In the case of Hinode/XRT, Figure [T71 shows the cor- 
relation between the DEM-weighted mean temperature [Tmean] and filter-ratio 
temperature [T rat i ] for four DEMs (regions) summarized in Tabled The T rat i 
is calculated for all suitable filter pairs which are investigated in Section I5T21 and 
summarized in Figure 1141 The values of T ra tio derived with different filter pairs 
are different, i.e. the thinner and thicker filter pairs give the lower and higher 
T r atio, respectively (see the "v" marks in Figures fT51 and [T4"|) . Additionally, T ra tio 
is higher than T m e a n in the lower temperature region (T < 1 MK: coronal hole), 
and T ra ti is lower than Tmean in the higher temperature region (T > 10 MK: 
flare). This is caused by the bias of XRT temperature response (Figure [TU]) 
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which has a sensitivity to mainly 1-10 MK plasma. However, the difference 
between T mean and T rat i is less than 50% (within two dotted lines in Figure [TT]) 
for all filter pairs. Hence, we claim that the filter-ratio temperature with XRT 
(the broad-band instruments) is useful for the quantitative analysis, e.g., coronal 
energetics. Additionally, the filter-ratio method has the following advantages: 

i) the filter-ratio method applied to XRT data can estimate a coronal tem- 
perature over a wide field of view at once, in contrast to a spectrometer. 

it] the filter-ratio method can investigate rapid temporal evolution of coronal 
temperature with a time scale well below one minute, since XRT can quickly 
take a data set with a pair of filters. 

Hence, the temperature diagnostics obtained with the filter-ratio method and 
the advanced DEM analysis with data sets taken with multiple filters of the 
XRT or EUV Imaging Spectrometer (EIS) onboard Hinode are complementary 
to each other. 



6. Summary 

We calibrated the effective area and coronal temperature response of the XRT 
using ground-based test data and also on-orbit data in X-rays and visible light, 
with the effect of contamination taken into account to the best of our knowledge. 
The calibrated thicknesses of pre-filter and FPAFs are summarized in Table [TJ 
The time- varying thicknesses of the contaminant accumulated on the CCD and 
FPAFs are shown in Figures [5] and [JJ respectively. 

Here, we examine how much our calibration results improve the coronal- 
tcmpcraturc-diagnostic capability of XRT. Figure [15] shows a comparison be- 
tween filter-ratio temperatures with our calibrated results and with the instru- 



mental parameters before calibration (Golub et al., 2007) for some suitable filter 



pairs (see Figure [H)) . This plot is created as following: We derive the expected 
intensity ratios with our calibrated parameters at each temperature [T ca \] in the 
useful temperature range of each filter pair (see Figure [H]) . Next, we derive the 
temperature [Tb c f] with the above-derived ratios and instrumental parameters 
before calibration. The vertical axis shows the fraction of the difference in the 
temperature estimate defined as (Tb f — T ca \) /T C3 \. Note that this plot varies with 
the accumulation of contaminants on the CCD and FPAFs as a function of time. 
This plot shows the case when one month passed after the CCD bakcout, i.e. 
800 A of contaminant was deposited on the CCD. While some filters happen to 
provide consistent temperature estimates in some temperature ranges, for most 
filters there are large deviations as illustrated in Figure [T51 Hence, we claim 
that our calibration result is indispensable for accurately estimating coronal 
temperatures with XRT. 

On the basis of our calibration results, we reviewed the coronal-temperature- 
diagnostic capability of XRT with the filter-ratio method. XRT has a wide variety 
of suitable filter pairs to estimate coronal temperatures with the filter-ratio 
method. For reference, we briefly mention the multi-temperature diagnostics. 
Using simultaneous data from more than two filters allows one to obtain more 
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Figure 18. Comparison between filter-ratio temperatures with the calibrated results in this 
article [T ca i] and with the instrumental parameters before calibration pb e f]. The horizontal 
and vertical axes of this plot show T ca i and the fraction of the difference in the temperature 
estimate defined as (7b c f — T ca i)/T ca i, respectively. The dotted line indicates the position where 
these temperatures are the same. 



information on the temperature distribution, and to resolve ambiguity in the 
results of the filter-ratio method (i.e., it increases the ability to pass Condi- 
tion (C) in Section [HIl]) . With many XRT filters, at least more than four filters 
( |Golub et al., 2007| , one can try to derive the differential emission measure 
(DEM). In this situation, the calculation of DEM with XRT data including 
the error estimate itself is the topic of interest (Schmcl z et ail I2009al [2009b; 
I2010p . The XRT is a powerful instrument not only to observe the morphology 
of coronal structures in detail but also to derive coronal physical quantities. 

6.1. Softwares for the coronal temperature diagnostics with XRT 

Lastly, we itemize the software used for the XRT calibrations and the coronal 
temperature diagnostics described in this article. These routines are distributed 
in the hinode/xrt directory tree of Solar Software ( jFreeland and Handy, 1998] ). 

xrt_ti me2conta m . pro 

gives the total thickness of contaminant on the FPAFs and the CCD as a 

function of time (Figures [5] and [7]) . 
xrt_eff_area.pro 

gives the XRT effective area (Figure O . 
xrtJIux.pro 

gives the XRT temperature response (Figure [TU]). 
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xrt_cvfact.pro 

gives the conversion factor [K^> (T)] from observed DN to incident photon 
number (Equation (fTU| and the top panel of Figures IT2"j) . With the /error 
keyword, the conversion factor [K^ (T)] from observed DN to the square 
of DN error is given (Equation (jT2J) and the bottom panel of Figures IT2"j) . 
xrt_teem.pro 

returns the estimated coronal temperature and emission measure with the 
filter-ratio method (Figure [TTj) . 

The XRT team have provided more software than the above. For example, 

ma ke_xrt_wa ve_resp . pro 

gives the XRT effective areas, and provides option for user-specified spectral 
emission models. 
make_xrt_temp_resp.pro 

gives the XRT temperature responses, and provides option for user-specified 
spectral emission models. 

The above two softwares are useful to examine the effect of uncertainties in the 
atomic physics, i.e. difference in spectral emission models, on the temperature 
analysis. For the DEM analysis, 

xrt_dem_iterative2.pro 

gives DEM solutions for multi-filter datasets, and provides option for Monte 
Carlo-style analysis of uncertainties. 

The details how to analyze the XRT data with those softwares are described in 



the XRT analysis guide (http://xrt.cfa.harvard.edu/resources/documents/XAG/XAG.pdf ).| 



Appendix 

A. Ground-Based Calibration 

Measurements of the X-ray transmission of focal-plane analysis filters (FPAFs) 
were made with five characteristic X-ray lines as part of the end-to-end through- 
put test with the flight XRT telescope (except for the flight pre-filters and 
sun-shield plate in front of the telescope) performed at the X-ray Calibration 
Facility (XRCF) at NASA/Marshall Space Flight Center in June 2005. Spec- 
tra of the X-ray lines used for measurement are shown in Figure [19] (see also 
Appendix lEl for details). 

A.l. Check of XRCF X-Ray Spectra 

In the end-to-end test, XRT pre-filters were not installed. Therefore, the effective 
area as measured at the XRCF can be written as 

A;ff@XRCF = A x R M1 x R M2 x Tfpafi x 7fpaf2 x QEccd- (17) 
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Figure 19. Spectra of characteristic X-ray lines and continua emitted by XRCF. These are the 
incident spectra into XRT. The solid lines are the spectra of continuum bremsstrahlung. The 
crosses (x) show the wavelength (horizontal axis) and photon- number flux in the unit of [cm -2 
sec -1 ] (vertical axis) emitted by the characteristic X-ray lines. If we indicate the characteristic 
X-ray lines as spectra in units of [cm -2 sec -1 A -1 ], the spectra can be indicated with the 
dotted lines, whose width and hight are 1 A and the same value as photon-number flux in the 
unit of [cm -2 sec -1 ], respectively. The gray areas show the error bar of the estimated spectra. 
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Figure 20. Left panel shows the residual between the DN measured with XRT CCD 
(-D./V mC as) and DN expected with XRCF spectra (DN cxp ) with one second exposure in the case 
where FW1 and FW2 are set at the open position. Black marks (— ) indicate the differences 
defined as -Rdiff = {DN cxp — D jV mcas ) / D Nm Ca , B . The top and bottom values near the black 
marks (— ) shows the value of DN cxp and DN mcSjB , respectively. The gray bar is the 1<t error 
bar. Right panel is the best correlative case, where the systematic error of 1.16 is taken into 
account. 



With FW1 and FW2 set to the open positions, this can be simplified as 

^off@XRCF = ^ X ^Ml x RM2 X QEqCD- (18) 

With the manufacturer-supplied information on the geometric aperture area 
[A], reflectivities of primary and secondary mirrors [Rmi an d R-M2.], and the 
Quantum Efficiency [QEqcd] of the CCD, we can check whether the XRCF X- 
ray spectra (Figure [H)]) were well estimated, by comparing signal from the XRT 
CCD (in terms of data number (DN) which represents the number of electrons 
generated) at the open/open position and that expected from energy spectra of 
the X-ray beams at XRCF measured with a flow proportional counter (FPC) 
(see Appendix lD.fi for details). 

The left panel of Figure [50] shows the result of the comparison. The verti- 
cal axis shows fractional difference i?diff = (DN cxp — DN mcas ) /DN meas , where 
DN meas and DN cxp are the measured DN with the XRT CCD and the expected 
DN from XRCF spectra. If DN meas and DN exp are consistent, i?diff ~ 0. The 
la error is shown by gray bars with observation error of CCD (e.g., photon noise 
and fluctuation of dark level) and the estimate error of XRCF spectrum taken 
into account as the source of errors. We found that i?diff of each X-ray line is 
almost constant at —0.14 as shown in left panel of Figure EOl This suggests that 
the spectral shape of XRCF beam is well estimated, but that the intensity is 
underestimated. We found that "DN exp x 1.16" gives the best fit to the data 
with Xdof=4 = 0-^9 as shown in the right panel of Figure [20j where DOF 
stands for the degree of freedom. We consider that this coefficient of 1.16 is 
the systematic error caused by the ambiguous aperture area of FPC and slight 
difference in the locations of XRT and FPC (see Figure l4"Tj) . because we do not 
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know the exact values of them and we apply the approximate values provided 
from XRCF. Hereafter, we derive DN cxp correcting this systematic error. The 
fact that DN mcas and DN exp show good match across the entire wavelength 
range in the right panel of Figure gives support to the reliability of the 
calibration result on the thicknesses of FPAFs with the XRCF X-ray spectra. 

A. 2. Calibration of Focal-Plane Analysis Filters 

Among the materials used for FPAFs, aluminum, titanium and beryllium can 
form thin oxidization layers on their surfaces. Hence, the effect of oxidization 
should correctly be taken into account when calibrating FPAFs. Oxidization of 
filters decreases transmission of FPAFs, especially at longer wavelengths, because 
the added oxygen absorbs X-rays. Considering oxidization, the transmission of 
FPAFs (7fpafi and 7fpaf2 for FW1 and FW2, respectively) arc expressed as 

7FPAF = 7purc X Tox x 7^ U pp , (19) 

where Tpurej Tox-, and 7^ U pp are the transmission of pure metal (unoxidized 
metal), oxidized metal, and support material, respectively. The transmission 
[T] of the material is determined by the material and its thickness [d] as 



r =" xp lruwJ- (20) 

where Z at t (A) is the attenuation length for the material at the X-ray wavelength 
A. The attenuation length of respective atoms (Z = 1-92, where Z is the atomic 
number) at an energy range of 50 - 30000 eV is derived by |Henke, Gull ikson, and Davis|| 

(|1993p and can be obtained from the "X-Ray Interactions With Matter" (http: / /hen ke.lbl.gov/optical_constants 
Using Equations (JT7J), C_g|), (H3), CP and also Equation (J31J) in Appendix [TJ7TT 
DN cxp at an exposure time of t = 1 sec is written as 



D7V exp = / (P (A) x ^ n XRCF (A) x exp ( ) x T ox (A) 



U (A) 



x7^ up p (A) x — ) dA, (21) 



where rf puro is the pure metal thickness of FPAFs, and P is the photon-number 
spectrum of the X-ray beam at XRCF (sec Figure IT9"1) . C gives the conversion 
from number of photons at wavelength A to signal DN generated on the CCD, 
and is expressed as C — t x /ic/e/3.65/G, where G is the system gain of the CCD 
camera, which is 57.5 [e DN" 1 ] for XRT ( |Kano et ai, 20080 . In performing the 
fitting, we set the support thickness d SU pp; namely 7^ U pp, as a fixed parameter. 
The reason for this is described below. For T ox , we adopt the oxidized- metal 
thickness derived with the on-orbit data, because after the end-to-end test at 
the XRCF and until the launch of Hinode, FPAFs were continuously purged 
with dry nitrogen to prevent additional oxidization. 

On the basis of the analysis in Appendix IB. 3. 11 the oxidized- metal thickness 
for the polyimide side can be regarded as A while that for the metal side (ex- 
posed to the ambient atmosphere) can be consistently set to 75 A, irrespective of 
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Table 4. Materials used in our calibration 



material 


molecular weight 


density 


Al 


26.98 


2.699 g cm" 3 


c 


12.01 


2.2 g cm" 3 


Be 


9.01 


1.848 g cm" 3 


Ti 


47.87 


4.54 g cm" 3 


A1 2 3 


101.96 


3.97 g cm" 3 


Ti0 2 


79.87 


4.26 g cm" 3 


BeO 


25.01 


3.01 g cm" 3 


polyimide (C22H10N2O5) 




1.43 g cm" 3 



the metal material. We note that our 75 A oxidized-metal thickness is consistent 
with the measurement of a well-oxidized aluminum filter bv lPowell et aZ.I (|1990|> . 
Hence, in this calibration, the fitting parameter is the metal thickness [d me tai] of 
each FPAF only. 

The results of fitting for each FPAF are shown in Figure [5T] and Figure 
For all filters, the fitted value [i?diff] is determined within the la error. The 
calibrated thicknesses of the FPAFs are summarized in Table \T\ 

In this fitting analysis, there are three notes: 

i) The thickness of polyimide film, which is used in Al-poly, C-poly, and 
Ti-poly filters as supports of these thin metal filters. According to the 
filter manufacturer (Luxel), the thickness of polyimide was well controlled. 
Therefore we fix the polyimide thickness to the values provided by the 
manufacturer (see Tabic [I]). 

ii) The fraction of the open area of the support mesh for the Al-mcsh filter. 
On the basis of the microscope measurement performed by the manufac- 
turer, we identified that the geometrical open area of mesh is 77%. In the 
calibration of Al-mcsh filter thickness, we adopted the value of 77% with 
scattering from the mesh pattern and the X-ray shape on it taken account. 
The detail of annulus transmission is described in Appendix [F] 

Hi) The last note is about the calibration of thick filters, namely mcd-Bc, 
thick- Al, and thick-Be filters. Since low-energy X-rays do not have enough 
transmission through the thick filters, we were able to use one or two X-ray 
lines to calibrate these thick filters. Thus we consider that the calibration of 
the thick filters is less certain. However, according to Equation ([20]) . even if 
the calibrated thicknesses of the thick filters have some error, the effect on 
their transmission is small. Hence we accept the thicknesses derived with 
XRCF data for now. To obtain further knowledge of the thicknesses of the 
thick filters, we need to take enough data sets with intense X-ray sources, 
i.e. active regions and flares, with thick filters. Calibration with on-orbit 
observation data is a project for future work. 

In Tablc[TJ wc sec that, for each metal filter, the difference between our thick- 
ness and that from the manufacturer is almost within the uncertainties reported 
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Figure 21. Calibrated thicknesses of focal-plane analysis filters mounted on filter wheel 1. 
The thicknesses described in the left, middle, and right side of "/" show the estimated pure 
metal, oxidized metal, and support thicknesses, respectively. Black marks (— ) indicate the 



differences defined as Rdiff = (DN C 



-DN n 



s)/DN u 



where DN„ 



is the DN measured 



with XRT CCD and D7V cxp is the DN expected with XRCF spectra. The top and bottom values 
near the black marks (— ) show the value of DN cxp and D./V mcas with one second exposure, 
respectively. The gray bar is the la error bar. 
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Figure 22. Calibrated thicknesses of focal-plane analysis filters mounted on filter wheel 2. 
The meaning of this figure is same as Figure [2T1 



by the manufacturer except for the Al-poly (difference is 200 A) and C-poly 
(800 A). In the case of the Al-poly filter, the manufacturer found a discrepancy 
between their measured thickness and that inferred from their own transmission 
measurement in visible light, and noted that the Al-poly filter should be thicker 
than their reported value. This information is consistent with our result. In the 
case of the C-poly filter, the difference is significantly large. This could have 
been caused by the different method of measurement. As we mentioned, we 
derived the thickness with the transmission in X-rays and attenuation length. 
The attenuation length depends on the density of material. In our calibration, 
we used the typical density of material summarized in Table U with the typical 
density of graphite adopted for our analysis on C-poly filter. However, it often 
happens that the assumption of a typical density is different from the actual 
density for the filter, especially in the case of C-poly filter, because carbon is not 
a metal. According to the manufacturer of filters, they have estimated that the 
density of their carbon film is about 90% of typical value of graphite, namely 
about 1.9 g cm -3 , because the carbon is evaporated using an electron beam in the 
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Figure 23. Thermal environment of XRT. The red, green, and blue lines indicate the 
temperature at the mirror support plate (XRTD_TEMP_11), at the focal-plane shutter motor 
(XRTD.TEMP.22), and at the CCD (XRTE.CCD.TEMP), respectively, where XRTD.TEMP.il, 
XRTD.TEMP.22, and XRTE.CCD.TEMP are names in the Hinode telemetry. 



process of making thin carbon filter. Considering this difference in density, the X- 
ray transmission of the C-poly filter estimated with our calibrated thickness and 
typical density is consistent with that calculated with manufacturer's thicknesses 
and w 90% of typical density because our calibrated thickness (5190 A, see 
Tabled]) is about 90% of manufacturer's thicknesses (6038 A). 



B. On-Orbit Calibration 

B.l. Thermal Distribution inside XRT 

Figure [23] shows the temporal evolution of the XRT thermal environment. The 
red, green, and blue lines indicate the temperature at the mirror support plate, 
at the focal-plane shutter motor, and at the CCD, respectively. From thermal 
viewpoint, there are three distinct episodes in XRT operation which arc shown as 
Phases 1, 2, and 3 in Figure [231 Phase 1 is the period from the launch of Hinode 
to first light, where the CCD bakeout heater was kept on to avoid contaminants 
accumulating on the CCD. When the CCD bakeout heater is on, we can say 
that no contaminants accumulate on the CCD, since the temperature of CCD 
becomes about 30°C, which is higher than its surroundings. Phases 2 and 3 are 
the periods of normal operation of XRT without and with enabling (turning on) 
the operational heater for the telescope tube, respectively. With the operational 
heater enabled, the rear end of the XRT telescope tube (in front of the CCD 
camera when viewed from the Sun) is warmed up to about +20°C. In these 
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periods, the CCD was cooled down to a temperature of about — 60°C to reduce 
the dark noise. Since the CCD is the coolest component in XRT, it is possible that 
the contaminant selectively accumulates on the CCD. The operational heater 
was originally kept off in Phases 1 and 2 due to technical reasons, and was then 
enabled from 18 June 2007. The period after this date is labeled as Phase 3. 

Let us now examine the relationship between the thermal environment of XRT 
and observed decrease in quiet-Sun X-ray intensity. The bottom panel of Fig- 
ure [5] shows temporal evolution of the X-ray intensity ratio normalized with the 
intensity observed with the Ti-poly filter to remove possible variation in emission 
measure for the quiet Sun (see Equations (|40| or (|4Tj) in Appendix ID.2|) . In the 
subsequent analysis, we made an assumption that the profile of the differential 
emission measure in the quiet Sun is almost stable, and the variation in the 
observed quiet-Sun intensity ratio can be attributed solely to the accumulated 
contaminant. During Phases 2 and 3, the CCD temperature was almost constant 
at around — 60°C. The decrease in the intensity ratio was continuously observed 
in these phases, indicating that the contaminant kept accumulating on the CCD 
in both Phases 2 and 3 (except the bakeout periods). The X-ray intensity ratio 
is recovered to almost the same level by each CCD bakeout as seen in Phase 3 
in the bottom panel of Figure [5J This implies two things: 

i) Each bakeout reduces the contaminant on the CCD to the same thickness 
(except for the "spots" discussed below). 

ii) Thickness of contaminant accumulated on each FPAF remains unchanged 
during Phase 3. 

The result i) is supported by the behavior of G-band intensity where the G- 
band intensity returned (reduced) to almost the same level by the CCD bakeouts 
(see Figure [25]). As will be mentioned in Appendix IB. 3.3[ we found that on/off 
operations of the operational heater when the temperature around the FPAFs 
is cooler than 5°C cause accumulation of contaminant on the FPAFs. Hence, 
we suspect that the accumulation of contaminant on the FPAFs took place in 
Phase 2 while not in Phases 1 and 3. 

Meanwhile, the first and third bakeouts created lumps of contaminant which 
are distributed over the whole CCD as spot-like patterns (see Figure HH in 
Appendix [Gj. These spots of contaminant cannot be removed by the subsequent 
bakeouts. We consider that excess amount of contaminant accumulated on CCD 
before the first and third bakeouts caused these spots. The detail of the spots is 
described in Appendix [G] In this regard, we consider that the contamination 
on the CCD comprises two seemingly different components: One is laminar 
contaminant that accumulates across the CCD as time goes on and can be well 
removed by a CCD bakeout when the accumulation thickness is small enough. 
The other is "condensed" contaminant forming spot-like patterns over the CCD. 
This condensation of contaminants occurs in association with a CCD bakeout 
when the accumulation thickness is large (threshold thickness is not known) and 
can not be removed by any further CCD bakeouts. In this article, we charac- 
terize only the laminar contaminant for contamination on the CCD. Detailed 
characterization of the spot-shaped contamination is a subject for future study. 
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Figure 24. G-band intensity enhancement observed after the second CCD bakeout from 
512 X 512 pixels at CCD center. The first G-band data just after the second bakeout has been 
taken as the reference for intensity enhancement. The gray areas in the left and right sides of 
this plot are the second and third CCD bakeout periods, respectively. 

We consider that no contamination accumulated on the pre-filter, because 
the temperature of the front end of the telescope (around mirrors) has been the 
warmest portion inside the telescope with the temperature at around or even 
exceeding 40° C in all phases. 

B.2. Visible Light Analysis - Contamination on the CCD 
B.2.1. Method 

We explain our method to measure the thickness of contaminant accumulated on 
the CCD using visible light (G-band) data. Figure |2~41 shows temporal evolution 
of G-band intensity after the second CCD bakeout on 3 September 2007. This 
plot is made from average G-band intensity over 512 x 512 pixels area at the 
center of the CCD for full-Sun (synoptic) images where the Sun center is located 
at the same position of the CCD (near the center of the CCD). Notable features 
in this plot are i) the intensity shows enhancement (non- negative increase) from 
the first G-band data just after the second bakeout, and ii) the intensity shows 
periodic oscillation. The period of oscillation is about 53.3 days. The maximum 
intensity enhancement is about 40% of the reference intensity. This amount of 
enhancement cannot be explained by the fluctuation of solar intensity which is 
less than 1% in visible light during the 11 year solar cycle. 

Our interpretation of the oscillation is that is caused by interference within 
a contaminant layer (anti-reflection effect). Figure [231 illustrates anti-reflection 
with a contaminant layer on the CCD. Both the direct ray (white arrow) and 
the reflected rays (black arrow) enter the CCD. Assuming three layers, vacuum, 
contaminant, and CCD, transmission of the visible light [7] into the CCD is 
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Figure 25. Anti-reflection effect by the contaminant accumulated on CCD. 
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(see Appendix [H] for the derivation of this equation) , where I is the intensity 
in vacuum, J t transmitted intensity into CCD, n refractive index of each layer, 
&2 wavenumber of G-band in the contaminant layer, and d thickness of the 
contaminant layer. The enhancement E of the G-band intensity is then given as 



T-T d=0 
Td=o 



sin 2 k^d x 



rc| (ni + n 3 )' 



(23) 



This equation clearly shows that the increasing thickness [d] of contamina- 
tion causes periodic intensity oscillation. Thickness of contaminant for a local 
minimum in the intensity enhancement is given by 



TO7T 

17' 



(24) 



where m is a positive integer (m = 1, 2, 3, • • •). The period of the intensity 
oscillation can be related to the contamination thickness as 



AT _ ]j_ _ ^2 _ 



2n 2 



(25) 



where ri2 is a refractive index of the contaminant and Ao is the wavelength of 
the incident visible light in vacuum (4300 A in G-band). This means that if we 
know the refractive index [772] of contaminant material, then we can measure 
the thickness of contaminant using the measured period of G-band intensity 
oscillation. 

In the above, we considered only contaminants on the CCD. But it is pos- 
sible that the contaminants also accumulated on the G-band glass filter. (Note 
this is indeed the case for the metal filters; see Appendix IB.3.1h Thus, the 
possible effect of contaminants accumulated on the G-band glass filter should 
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Figure 26. Expected amplitude (maximum value) of the intensity enhancement by the 
anti-reflection effect (a in Equation I I26I I) is shown by contours. The observed amplitude of the 
intensity enhancement (0.444; see Figure I24I) is indicated by a thick contour. The refractive 
indexes of CCD (713 ks 4.9) and G-band filter (713 ss 1.6) are shown in orange and green, 
respectively. 

be considered. As will be detailed in Appendix IB.3.11 we found that the con- 
taminant is a long-chain organic compound which consists of carbon, hydrogen, 
and oxygen, and that its refractive index is close to 1.5. With this information, 
let us next consider the anti-reflection effect with the contaminant accumulated 
on the G-band filter. Figure [55] shows expected amplitude (maximum value) of 
intensity enhancement caused by the anti- reflection effect with Equation ([2~3l) . 
The material (refractive index) of G-band filter, CCD and contaminant are glass 
(77,3 w 1.6), silicon (713 w 4.9) and contaminant (77,2 ~ 1-5), respectively. 
As seen in the green area of Figure 1261 intensity enhancement caused by (pos- 
sible) contaminant on the G-band filter is close to zero due to the proximity of 

77.2 ~ 1.5 and 77,3 « 1.6. Meanwhile, the contaminant on the CCD can enhance 
the G-band intensity enough (more than 0.4) for the observed enhancement 
shown in Figure [Ml Hence, even if the thickness of contaminant accumulated 
on the G-band filter changes as a function of time, we can ignore the G-band 
intensity enhancement caused by the contaminant on the G-band filter. 

By assuming that the refractive index of the contaminant is 1.5 and from 
Equation ([25]) . one period of intensity oscillation corresponds to 1433 A thickness 
of contaminant as an average over 512 x 512 pixels at the center of CCD. The 
accumulation rate of the contaminant there is estimated to be 9800 A year -1 , 
because the period of the G-band intensity oscillation at the CCD center is about 

53.3 days. 

B.2.2. Model 

In order to avoid excess accumulation of contaminant on the CCD, we started to 
perform regular CCD bakeout every three -four weeks since the third bakeout. 
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Figure 27. Relation between intensity enhancement in G-band and thickness of contami- 
nant on CCD. The dotted curves show the exponentially decaying amplitude of the intensity 
enhancement oscillation. 



This means that, for the period after then, there is not much time to see the 
intensity oscillation in the G-band profile between two adjacent bakeouts. Hence, 
a quantitative estimate on the contaminant thickness directly from the G-band 
intensity turned out to be not possible. Instead, for the period after the third 
bakcout, we developed an empirical model of the G-band intensity oscillation 
using the data between the second and third CCD bakeouts. On the basis of this 
model, we estimated time-dependent thickness of the contaminant on the CCD 
after the third CCD bakeout in addition to the period between the second and 
third CCD bakeouts. 

In the profile of the observed intensity enhancement (Figure [23]), the periodic 
oscillation and the decay of amplitude are seen. The empirical model incorporates 
these two properties. Modifying Equation (|23[) . the empirical model, which re- 
lates intensity enhancement in G-band (i? m odei) to the thickness of contaminant 
[d] , is expressed by the following equation whose profile is shown in Figure 1271 

x exp(-6 x d) + ^, (26) 

where a and b are the amplitude (maximum value) and decay coefficient of 
intensity enhancement, respectively. The inside of { } in Equation (|26[) is the 
rewriting of Equation (|23[) with 

g = ( n i - n D ( n i - n i) (27) 

1 + a n \ (n x + n 3 ) 2 

The best-fit values to the observed G-band intensity enhancement (see Figurc l24l) 
are a = 0.444 and b = 8.23 x 1CT 5 . 

In this empirical model, there are two notes: 

i) With the observed amplitude (maximum value) of intensity enhancement 
(a w 0.444) and the refractive index of CCD (713 w 4.9), the refractive 
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Figure 28. Temporal evolution of G-band intensity monitored with 512 X 512 pixels at CCD 
center. 
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Figure 29. G-band data (left panel) and estimated thickness of contaminant on the CCD 
(right panel) after the sixth bakeout. 



index of contaminant is calculated as « 1.4 (see Figure l2"fj]l . This 
value is certainly different from our assumption that the refractive index of 
contaminant is 1.5. 

ii) In the observed intensity profile (Figure IIM]) . there is slight decrease in the 
amplitude, which is incorporated as an exponential decay in the model (b 
in Equation ((26)) ) and as indicated by dotted lines in Figure [27] 

As will be discussed in Appendix IK.21 we conclude that they do not significantly 
affect the results of XRT calibration. 

We continuously monitor the temporal evolution of G-band intensity (Fig- 
ure 128]) . As seen in Figure [28l only a portion of an entire cycle of the intensity 
enhancement, namely initial rise of the intensity, is available for G-band data 
for each bakeout after 1 February 2008 (after the third bakeout). We apply 



SOLA: Narukage_etal_R3.tex; 15 November 2010; 1:16; p. 43 



N. Narukage et al. 



the empirical model given by Equation ([2l)]) to that period. The left panel of 
Figure [25] is the observed G-band intensity enhancement after the sixth bakeout, 
and the right panel shows the thickness of contaminant on the CCD derived 
from the model. The orange line in the right panel shows the linear least-square 
fit to the thickness, which yields an accumulation rate of contaminant on the 
CCD of about 613 A month -1 , where 30 days are taken as one month unless 
otherwise noted. In this article, we assume a constant accumulation rate for the 
contaminant between two successive bakeouts. 

We note that the X-ray intensity ratio and G-band intensity are recovered to 
almost the same level by each CCD bakeout as seen in Phase 3 in the bottom 
panel of Figures \8\ and [28l respectively. This suggests that after each bakeout 
the contaminant reduces to the same thickness. Furthermore, the recovered G- 
band intensity in Phase 3 is consistent with the G-band intensity taken just 
after the end of Phase 1 (until which we expect no contaminant on the CCD due 
to continuous CCD bakeout since launch). Hence, we conclude that this "same 
thickness" is actually zero. 

Table [7] in Appendix [Cl summarizes the thus-derived accumulation rate of the 
contaminant between each bakeout, and Figure [5] indicates temporal evolution of 
contaminant thickness on the CCD. The similar plots to Figure |2T)1 after all bake- 



outs arc distributed with Solar Software (SSW) as SSW_DIR/hinode/xrt/idl/response/contam/xrt_contam_on_cc( 

where SSW_DIR indicates SSW directory in your environment. This file and 
database of contaminant thickness will be regularly updated as the CCD is 
baked out repeatedly. 



B.2.3. Spatial distribution of the contaminant across the CCD 

To monitor the accumulation rate of contaminant on the CCD with G-band 
data, we have used the 512 x 512 pixels at the center of the CCD. On the 
other hand, when looking into fine structure of contaminant accumulation for 
the entire imaging area of the CCD, we note that there is certain spatial distri- 
bution of accumulation rate of contaminant across the CCD (see the left panel 
of Figure 130]) . The accumulation rate decreases towards the edges of the CCD 
from its center. The difference of accumulation rate between the CCD center 
and outermost area in the left panel of Figure IBTil is about 10%. 

By building a numerical contamination model in which geometry and tem- 
perature of structures around the CCD have been incorporated, [Urayama et al.\ 
(2008) studied the expected spatial distribution of contaminants accumulated 
on the CCD (see the right panel of Figure 130]) . The resultant distribution was 
quite similar to the observed one (see Figure I30[) when the contaminants come 
from the forward direction of the CCD, while not when the contaminants come 
from the side of the CCD where the pre-amplifier unit for the CCD is located. 
The similarity is seen not only in the distribution of the contaminant but also in 
the relative difference in the accumulation rate of about 10% between the CCD 
center and the outermost area in the observed data (left panel of Figure [30]) • 
Since the edge of the CCD is located outside the solar disk (see the left panel 
of Figure I30[) with synoptic G-band exposure, we can not measure the accumu- 
lation rate at the edges of the CCD. However, on the basis of the comparison 
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Figure 30. Observed and simulated spatial distribution of the contaminant across the CCD. 
The left panel is the observed accumulation rate map with G-band intensity oscillation. The 
background is a full-Sun (synoptic) image in G-band. The image size is the full CCD area 
(2048 X 2048 pixels). Each colored box corresponds to 128 X 128 pixels on the CCD. The right 
panel shows the simulated accumulation-rate map under the condition where the sources of 
contaminant located in the forward direction from the CCD. The simulated area in the right 
panel is also full CCD area. 



between the results of observation and numerical simulation, we conclude that 
the difference of accumulation rate between the CCD center and its edges is 
about 20%. 



B.3. Analysis of Contaminant with X-Ray Data 



Next, we examine whether the decrease in quiet-Sun X-ray intensity shown in 
Figure [5] is consistent with the effect of the contaminant accumulated on the 
CCD whose thickness derived from G-band data. For this analysis, we used data 
sets in the period between the second and third bakeouts (in Phase 3), because 
of availability of data sets both in G-band and in X-rays. 

Plus signs in Figure [3T] indicate temporal evolution of quiet-Sun intensity 
ratios after the second CCD bakeout made from sets of simultaneously taken four 
filter images: Al-mesh, Al-poly, C-poly, and Ti-poly. In order to remove possible 
variation of emission measure in different quiet-Sun regions, each intensity was 
normalized by the geometric mean of intensities observed with the four filters, 

-mesh X -^Al-poly X -^C-poly X -^Ti-poly) ^ ; (28) 

where /filter is the X-ray intensity with each filter. Data points in this plot were 
obtained from full-Sun images showing no active region in order to remove the 
effect of scattered light from active regions. Data gaps in this plot are due to 
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Figure 31. Observed and expected X-ray intensity ratio after the second CCD bakeout in 
case of no contaminant on focal-plane analysis filters. 



the appearance of active regions. In the plot, a significant decrease in intensity 
ratio with the Al-mesh filter and a slight increase in those with the other three 
filters arc clearly seen. We consider that these changes are caused by the effect 
of accumulating contamination, which absorbs X-rays. Note that the increase 
in intensity ratio with Al-poly, C-poly, and Ti-poly filters is due to the large 
decrease in intensity with Al-mesh filter (see the top panel of Figure [5]). 

Next, we calculate the expected intensity of the quiet Sun observed with the 
contaminant-accumulated CCD using the differential emission measure (DEM) 
at the quiet Sun. From Equation (fT]), the effective area of XRT ^ontam^ including 
the effect of contamination, is given by 

^contam _ x 7^ ontam (29) 

where 7^ ntam is the transmission of the contaminant. For the DEM of the quiet 
Sun, we adopted the DEM profile derived bv lBrooks and Warrenl (2006 ) with an 
assumption that the profile of the quiet-Sun DEM is almost constant in time. 
The procedure for deriving the expected intensity from the thus-assumed DEM 
is detailed in Appendix ID. 21 We will discuss the effect of the difference in the 
DEM model on our analysis in Appendix IK. 41 

By adopting the thickness of contaminant on the CCD (from the G-band data) 
for deriving % ntam, the combination of ^.contam an( j DEM profile gives the 
expected profiles of the intensity ratio, which are overlaid in Figure 1311 (The 
chemical composition of DEHP (Diethylhexyl phthalate: C24H38O4) is adopted 
as the contaminant on the CCD for making the profiles in the figure; see Ap- 
pendix [B.3.1l for details.) The figure clearly shows that the observed and expected 
X-ray intensity ratios are inconsistent with each other. Especially, in addition 
to the large deviation in Al-mesh filter, although the expected intensity ratio 
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Table 5. Candidates of contaminant material 
# material density refractive index 

1 diethylhexyl phthalate^ 0.986 g cm~ 3 n = 1.5^) 

(DEHP) C24H38O4 

2 tctramethyl tetraphenyl trisiloxane' c ' 1.07 g cm~ 3 n = 

(MPS) C28H 3 2 2 Si3 

3 polydimethylsiloxane( c ) 0.971 g cm~ 3 n = 

(PDMS) (C 2 H e OSi)„ 

( a ) DEHP is a plasticizer. 

W catalog of chemical reagent by Wako Pure Chemical Industries, Ltd 

< c ) MPS and PDMS are classified Siloxane. 

M IQsantowskil (fT982l > 

< e ) catalog of TORAY Silicone, Ltd 



with the Al-poly filter is larger than that for C-poly filter, the observed intensity 
ratios with these two filters arc contrary to expectation. In the following, we 
investigate the cause of this inconsistency. 

B.3.1. X-ray Analysis for Phase 3 

In order to solve the intensity controversy between Al-poly and C-poly filter, 
in this section, we characterize the contaminant accumulated on each FPAF. 
We note that contamination on each FPAF is the only possible candidate to 
account this controversy, because the other optical elements are common for 
the data taken with any of the X-ray filters. As discussed in Appendix IB. 11 it 
is expected that there was no additional accumulation of contaminant on the 
FPAFs in Phase 3. As the period shown in Figure [3T] belongs to Phase 3, the 
data shown in the figure is suitable for identifying thickness of contaminant on 
each FPAF. Moreover, using these data, we can also identify the properties of 
contaminant material, i.e. the chemical composition, density and refractive index 
of contaminant material, because the transmission of the contaminant [7^ ntam] 
in Equation ((29|) depends on the above properties as shown by Equation (l20l) . 
where the thickness [d] of contaminant is estimated with refractive index of 
contaminant material (see Equation ([^5]) ). and the attenuation length [Z a tt (A)] 
is determined by the chemical composition and density. 

Before proceeding to the above calibration, we first try to identify the possible 
material of the contamination. We picked three candidate materials, which arc 
listed in Table [5J These three candidates are widely used for satellites, and are 
well known as possible source of contamination: DEHP (Diethylhexyl phthalatc: 



SOLA: Narukage_etal_R3.tex; 15 November 2010; 1:16; p. 47 



N. Narukage et al. 



Table 6. Fitting parameter to observed X-ray data 



parameter 


examined range 


expected value 


best fit value 


contaminant 


DEHP, MPS or PDMS 




DEHP 


contam. on Al-mesh filter 


0-4000 A 




1200 A 


contam. on Al-poly filter 


0-4000 A 




2900 A 


contam. on C-poly filter 


0-4000 A 




500 A 


contam. on Ti-poly filter 


0-4000 A 




400 A 


metal thick, of pre-filter 


1000-2000 A 


1538 A< a > 


1550 A< b > 


oxide thickness' ) 


0-75 A 


< 75 A« 


75 A / o A( c ) 



( a ' This value is based on the certification sheet (sec Figure [3] and Table \Q. 

( b ' The expected metal thickness when it was fabricated. At fabrication, the metals had not 
oxidized at all. The on-orbit pure and oxidized- metal thicknesses are shown in Table [T] 

( c ) The oxide thickness of metal filters. We examine it for open side and for polyimidc side. 

W The oxide thickness of well-oxidized Aluminum filter is 150 A |Powell et a l., 1990] > . Hence, 
the expected oxide thickness for one side (open side) is < 75 A. 

( e ) Best fit value of "oxide thickness for open side" / "oxide thickness for polyimide side" . 



C24H38O4) is the representative of a long-chain organic compound which consists 
of carbon, hydrogen, and oxygen. MPS (tetramethyl tetraphenyl trisiloxane: 
C28H32O2S13) and PDMS (polydimethylsiloxane: (C2H 6 OSi)„) are long-chain 
organic compounds including silicon. Although these candidate of contaminant 
classified into two types, namely a long-chain organic compound without and 
with silicon, their density and refractive index are almost the same: 1 g cm -3 
and 1.5, respectively. We note that, though DEHP is not used in the XRT, 
if the actual contaminant material is a long-chain organic compound without 
silicon, the chemical composition, density, and refractive index would be similar 
to DEHP. Hence, in our analysis, we refer to DEHP as the representative of a 
long-chain organic compound without silicon. 

We obtained filter parameters that best fit to the observed intensity ratio. The 
parameters and their examined range for the fitting are summarized in Table [5] 
The examined range of contamination thicknesses accumulated on the four thin 
FPAFs is 0-4000 A, since contaminant thickness larger than 4000 A gave results 
that are inconsistent with the observed data. In this calibration, we cannot 
investigate the contamination thickness on thicker FPAFs, because the thicker 
FPAFs cannot observe the quiet Sun enough. The metal thickness of the pre-filter 
was examined in the range of 1000 - 2000 A, i.e. approximately covering the range 
of ± 500A from the thicknesses provided by the manufacturer. According to 
IPowell et al.\ (|1990|) , the oxide thickness of a well-oxidized stand-alone aluminum 
filter is 150 A in total for both sides, i.e. 75 A for each side. Hence, we examine 
oxide thickness for open side of a filter in a range of - 75 A under the assumption 
that any metal filters were oxidized with the same thickness. Some of the FPAFs 
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best-fit X-ray intensity ratio (DEHP assumed) 



Phase ] 
4 : 



3 - 



2 - 




1 1 

observed intensity ratio 

best-fit intensity ratio 
with CCD and filter 
contaminant (DEHP assumed) 



Al-mesh / 4 filters 
Al-poly / 4 filters 
C-poly / 4 filters 
Ti-poly / 4 filters 



Sep 



Oct Nov Dec Jan 

date (from 01-Sep-2007 to 31-Jan-2008) 



Feb 



Figure 32. Best-fit X-ray intensity ratio after the second CCD bakeout, where for this 
analysis we have taken the material to be similar to DEHP. 




Figure 33. Best-fit X-ray intensity ratio after the second CCD bakeout, where the materials 
of contaminant are assumed to be MPS (left panel) and PDMS (right panel). 



have a support made of polyimide. We also investigated oxide thickness for 
polyimidc side in the range of - 75 A under the above assumption. 

Figure [35] shows the estimated X-ray intensity ratio from the best-fit pa- 
rameters summarized in Table El including material for the contaminant, which 
was identified to be consistent with DEHP, namely the representative material 
of a long-chain organic compound without silicon. The area-weighted average 
thickness of the pre-filter was calibrated to be 1550 A, which is the expected 
metal thickness when it was fabricated. At fabrication, the metals had not 
oxidized at all. This value is very close to the value of 1538 A derived from 
the manufacturer-supplied information on thickness for each piece of the pre- 
filter (sec Figure [3]). Hence, considering the oxidization (using Equation (f5rJf of 
Appendix |T| , the on-orbit pure and oxidized- metal thicknesses of pre-filter are 
calibrated to be 1492 A and 75 A, respectively, as summarized in Table [TJ The 
calibrated thickness of contaminant on each FPAF turned out to be different. 
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observed X-ray intensity ratio after the second bakeout 



* 2 | relation between X-rays and contaminant 

7 1 r 1 — ' i ' — i ; 3000.9 7 E 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 




date (from 01-Sep-2007 to 31-Jan-2008) thickness of contaminant on CCD [A] 

Figure 34. left panel: Observed X-ray intensity ratio (black +) and measured thickness of 
contaminant on the CCD with G-band (magenta line) after the second bakeout. right panel: 
Relation between X-ray intensity ratio and thickness of contaminant on the CCD after the 
second bakeout. 



The thickness of contaminant on the Al-poly filter was identified to be 2900 A 
which was the thickest. We will see in Appendix IB.3.3I that this difference in 
thicknesses is consistent with the frequency of filter usage. 

Figurc[33]indicates difference in the fit among different candidates for the con- 
taminant material. The left panel shows the best- fit result with MPS, while right 
panel PDMS. Clearly, these materials are not able to account for the measured 
profiles, especially that for the Al-mcsh filter. Difference in chemical composi- 
tions among candidate materials for contamination most affect transmission at 
longer wavelength X-rays, hence intensity with the thinnest Al-mcsh filter. This 
results in the situation that the Al-mesh profiles show the largest sensitivity 
against to the different candidates. Furthermore, for MPS and PDMS, the metal 
thickness of the prc-filtcr which gives the best fit to the measurements turned 
out to be fOOO A and 1050 A, respectively. These thickness are unlikely because 
of significant deviation from the averaged pre-filter thickness measured by the 
manufacturer (1538 A). Hence, we can conclude that the contaminant material 
is a long-chain organic compound whose characteristics, namely chemical com- 
position, density, and refractive index, are similar to DEHP. We will discuss this 
identified contaminant material in Appendix IK. II 

B.3.2. X-Ray Analysis for Phase 2 - Contamination on the CCD 

In Appendix IB. 21 we established the method to monitor the contaminant on 
the CCD with G-band data. Using this method, we can measure the thickness 
of contaminant on the CCD for the Phase 3 period. However, for Phase 2, we 
cannot apply this method because there is not enough G-band data. Instead, 
we rely on the X-ray data for investigating accumulation of contaminant on the 
CCD. 

For this purpose, let us first derive the relationship between G-band and X- 
ray data for the period after the second bakeout (which took place in Phase 3) 
where both data are available. The measured thickness of contaminant on the 
CCD with G-band data is shown in magenta in the left panel of FigurcGHl In this 
panel, the observed X-ray intensity ratio between Al-mesh and Ti-poly filters is 
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observed data before the first bakeout 

Phase ZXZ^^^^^^^^^H^^^^^^^^^HZZl 



4.0 "I ' 1 ' 1 ' 1 ' 1 ' : 0.5 




2.0 U !Li I i I i I i I i 3 0.0 

Oct Dec Feb Apr Jun Aug 

date (from 22-Sep-2006 to 31-Jul-2007) 



Figure 35. Observed X-ray intensity ratio (black +) and intensity enhancement in G-band 
(magenta x) before the first bakeout. 

also shown in black color. Using these data, we obtain the relation to convert 
the X-ray intensity ratio to the thickness of contaminant on the CCD which is 
shown in the right panel. As demonstrated in Appendix IB. II the contaminant 
on FPAFs is constant in Phase 3. Note the necessary condition for applying the 
right-panel relation of Figure [3H (hereafter, "XSC relation", where XSC stands 
for X-ray-suggested CCD contamination) is that the contaminant thicknesses on 
the FPAFs are the same as those in Phase 3. 

Now, let us check whether the above necessary condition is satisfied in Phase 2 
too. Figure [35] shows the temporal evolution of X-ray intensity ratio (black +) 
and intensity enhancement in G-band (magenta x ) observed before the first 
CCD bakeout. Note the entire period of Phase 2 is contained in the plot. By 
applying the XSC relation to the observed X-ray intensity ratio, we estimated the 
thickness of contaminant on the CCD as shown by the black pluses in Figure [551 
Note the thickness was estimated by assuming the existence of contaminant 
on the FPAFs whose thicknesses are identical to those in Phase 3. Thickness 
of the contaminant on the CCD derived with the G-band data is also plotted 
in Figure [55] in magenta crosses. This profile has some trends. We classify the 
profile into A-E as shown in Figure 1551 In periods A, C, and E, the thickness of 
contaminant increases at a certain rate, but in B and D, the thickness is almost 
constant. As discussed earlier, the variation of G-band data is solely affected 
by the contaminant on the CCD, not by the contaminant on the G-band filter. 
Hence, using the G-band data, we can monitor the contaminant on the CCD 
independently. The magenta crosses in Figure [36] are the estimated thickness 
of contaminant on the CCD from G-band data. As clearly seen, in period B 
the thickness of contaminant on the CCD estimated from X-ray intensity ratio 
(black +) is consistent with that derived with G-band data (magenta x). This 
suggests that the above necessary condition of the XSC relation is valid. Now 
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Figure 36. Accumulation rate of contaminant on CCD before the first bakeout. This rate is 
derived from two kinds of data sets: X-ray intensity ratio between Al-poly and Ti-poly filters 
(black +), and G-band intensity (magenta x). The orange solid and dotted lines show the 
measured and expected thickness of contaminant on the CCD, respectively. The measured and 
expected accumulation rates are also shown with red and blue characters. The [A] — [E], [X], 
and [Y] show the classified periods from the profile of contaminant thickness (orange lines). 



that we know that the thicknesses of contaminant on the FPAFs in period B are 
the same as those in Phase 3 and that it is not likely the thickness would decrease 
in any certain period (note they did not decrease even with FPAF temperature 
~ + 20° C in Phase 3), we can conclude that the thicknesses of contaminant on 
the FPAFs are unchanged after (including) period B (after (including) the "late 
period" of Phase 2 in Figure [5]). This means that, after period B, the thickness 
of contaminant on the CCD can be measured from the XSC relation as shown 
by the black +. The solid orange lines after period B show a linear fit to the 
measured thickness (black +). 

For the interval after period B, there are episodes of rapid accumulation of 
contaminant on the CCD (periods C and E), while the rest do not show much 
accumulation (periods B and D). We note that the behavior of contaminant 
accumulation on the CCD is closely related to the status of the operational 
heater as follows: 

• From turning on the operational heater to about ten days after turning 
off, which we call "contamination period (CP)", the contaminant rapidly 
accumulated on the CCD at a rate of < 700 A month -1 (periods C and E 
but see comments below) . The periods when the operational heater was on 
are indicated by blue areas in Figure I3U1 

• From about ten days after turning off to the next turning on, which we 
call "small-contamination period (SCP)", the accumulation rate is small, 
namely < 100 A month -1 (periods B and D). 
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Hereafter, we call this relationship the "ODC relation" where ODC stands for 
operational heater driven contamination. This relationship suggests that the 
contamination is triggered by the operation of the operational heater. 

So far, we have seen that the XSC relation can be applied for the period from 
period B ("late period" of Phase 2 in Figure [5]). Here let us check if this relation 
still holds for the rest of Phase 2, namely period before (including) period A 
on Phase 2 ("early period" of Phase 2). From Figure [23 there are the following 
trends in the temperatures: 

i) Until March 2007, when the operational heater was on, the temperature 
around the FPAFs was ps + 20°C and at the CCD ps - 65°C. When 
the operational heater was off, the temperature around the FPAFs was 
ps + 5°C and at the CCD ps - 70°C. 

ii) From March 2007 to May 2007, the temperatures around the FPAFs and 
at the CCD increased a little due to the increased Earth albedo towards 
the beginning of the eclipse season. 

Hi) In the eclipse season (from May 2007), the temperatures around the FPAFs 
decreased due to the decreased time of the solar illumination, while the 
temperature at the CCD increased due to the increased Earth albedo. 

These changes in the temperatures were negligible for the accumulation rate of 
the contaminant on the CCD, since the measured accumulation rates were almost 
constant during periods B and D, which correspond to the above trends i) - 
ii) and trend Hi), respectively. Therefore, we expect the ODC relation is also 
applicable to the period before (including) period A. 

There are not sufficient X-ray nor G-band data for the interval between pe- 
riods X and A ( "early period" of Phase 2 in Figure [SJ to measure contaminant 
thickness on the CCD. Therefore, we rely on the ODC relation for this interval 
and check if the relation consistently connects the periods from A back to X. 
To apply the ODC relation, we assumed that the contaminant accumulated on 
the CCD at a rate of 700 A month- 1 for CP while 100 A month- 1 for SCP, as 
shown by dotted orange lines in Figure 131)1 In performing the extrapolation with 
the ODC relation, we started from period B and backwards toward period X. 
Note that the applied rates are somewhat larger than the measured rates in 
periods C (580 A month -1 ) and periods A and D (60 A month- 1 ). This is 
because we expect that the accumulation rate in early period in Phase 2 would 
be higher than the late period in Phase 2. Note here, that the difference in 
contaminant thicknesses at the beginning of period X is small (ps 200 A) even 
if we adopt the 580 A month -1 and 60 A month- 1 pair. For period Y, which is 
the interval between the end of Phase 1 (until which we expect no contaminant 
on the CCD due to continuous CCD bakeout since launch) and the beginning 
of period X, a constant accumulation rate of contaminant on the CCD was 
assumed. We note that the duration of period Y is determined to be ten days 
as shown in Figure 1361 since we confirmed that the X-ray intensity in the quiet 
Sun observed with Al-poly filter was almost constant for the period when more 
than ten days passed after the beginning of Phase 2. Note that there is no data 
available in the first ten days. The estimated accumulation rate in period Y is 
ss 1400 A month -1 . This high rate might be caused by the hot temperature 
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of the telescope in Phase 1. The temperature of the front-end portion of the 
telescope turned out to be significantly hotter than the temperature predicted 
before launch. We guess that in Phase 1 a great deal of contaminant was created 
and filled the inside of the telescope. This contaminant rapidly accumulated on 
the CCD after the CCD bake heater was turned off (in the beginning of Phase 2, 
i.e. that of period Y). 

Finally, let us check whether the expected contaminant thickness on the CCD 
is reliable. In periods X and Y, the expected thickness (dotted orange line) is 
consistent with the contaminant thickness measured with G-band data (magenta 
x). This result strongly supports our expectation. As discussed above, now 
we understand the accumulation of contaminant on CCD in all periods (see 
Figure IHJ). 

We note that in period A there is a discrepancy between the contaminant 
thickness on the CCD derived with the ODC extrapolation (dotted orange 
line) and the thickness derived with the XSC relation (black pluses +). This 
discrepancy will be discussed in the next section. 

B.3.3. X-Ray Analysis for Phase 2 - Contamination on the FPAFs 

From discussion in Appendix IB. II and Appendix IB.3.21 we know the following 
for the contamination on the FPAFs: 

• The contaminant did not accumulate in Phase 3. 

• In Phase 2, the contaminant did not accumulate after the beginning of 
period B. 

Hence we conclude that the FPAFs contaminant must have accumulated by the 
end of period A. 

Figure [37] shows temporal profiles of the temperature around the FPAFs 
and estimated thickness of contaminant on the CCD. In period A, there is a 
discrepancy between the thickness of contaminant on the CCD derived with 
the ODC relation (dotted orange line) and the thickness derived with the XSC 
relation (black pluses +). Given the successful ODC extrapolation described 
in the previous section, this discrepancy most likely implies that the thickness 
of contaminant on the FPAFs that we assumed (being the same as those in 
Phase 3) for applying XSC relation in period A is not correct. In other words, 
the contaminant accumulated on the FPAFs at least in period A. 

Next, we investigate why the FPAFs contamination took place in period A 
while not in period C. Both periods A and C are "contamination periods (CPs)" 
defined in Appendix IB.3.21 As we mentioned, in these CPs, the contaminant 
accumulated on the CCD. The difference between periods A and C is the tem- 
perature around the FPAFs as shown in Figure [27J The temperature around 
the FPAF in period A (about 4°C) was cooler than in period C (about 8°C). 
We classify the CPs into two according to the temperature around the FPAFs: 
For the periods when the temperature is lower than ps 5-8°C, we call "cool 
contamination period (CCP)", while "warm contamination period (WCP)" for 
the rest of the CP. The difference in temperature leads to the hypothesis that 
the contaminant generated by turning on the operational heater accumulated 
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Figure 37. Relation between the temperature around the focal-plane analysis filter (green 
line) and estimated thickness of contaminant on the CCD with X-ray intensity ratio (or- 
ange line). The contaminant would accumulate on FPAFs in the "cool contamination periods 
(CCPs)" shown in magenta. While, in the "warm contamination period (WCP)" shown in 
yellow, the contaminant would not accumulated. The blue dashed line indicates the boundary 
temperature between CCP and WCP (5°C). 



on the cool FPAFs in CCPs. On the other hand, the high temperature around 
the FPAFs in the WCPs prevents the contaminant from accumulating. If this 
hypothesis is correct, the contaminant accumulated on the FPAFs in the five 
CCPs shown in magenta in Figure [37J On the other hand, in the WCP shown 
in yellow, the contaminant did not accumulate on the FPAFs. Although the 
difference of temperature around the FPAFs is only an indirect evidence, there 
is not any reasonable timing for the FPAFs contamination other than the CCPs. 
On the basis of this result, we expect that the contaminant did not accumulate 
on the FPAFs in Phase 1, since the temperature around the FPAFs was higher 
than 10° C in that phase. 

As we identified that the CCPs are the periods when the FPAFs suffered from 
the contamination, next we examine the amount of the contaminant accumulated 
on each FPAF as a function of time for each CCP. Figure 1551 gives the frequency 
in use of FPAFs during the CCPs. The vertical axis shows the cumulative time 
for each filter located on the optical path. We expect that this cumulative time 
should have a positive, linear correlation with the thickness of contaminant 
for each FPAF. Hereafter, we call this relationship the "USF relation" where 
USF stands for FPAF-usage-suggestcd FPAF contamination. To estimate the 
thickness, we apply the following simple model: We set the accumulation rates 
of contaminant on FPAFs to be x, y, z, and w as shown in Figure [3U Here, x 
and z arc the accumulation rates on the filters located on the optical path (o.p.), 
while y and w are the rates on the filters located just next to the optical path 
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frequency in use of filters 
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date (from 22-Sep-2006 to 30-Jun-2007) 



Figure 38. Cumulative time of filter usage during the "cool contamination periods" colored 
by magenta in Figure [371 The final percentage of cumulative time for each filter is described 
in a parenthesis, where the Al-poly filter has the largest cumulative time. 




Figure 39. Modeled accumulation rate of contaminant on each focal-plane analysis filter, x 
and y are the rates of filters mounted on filter wheel 1, and z and w are rates of filters on filter 
wheel 2. 



(o.p. ± 1, in Figure 131?)) . We assume that no contaminant accumulated on the 
filters far from the optical path (o.p. ± 2 and o.p. + 3). Using these accumulation 
rates, we can express the thickness of contaminant accumulated on each filter d 
as follows: 

dFPAFl = X X * F PAF1 + y X (tppAF 1 ! + 4'pafi) i ( 30 ) 
^FPAF2 = Z X t F pAF2 + w X (*FPAF2 + *FPAF2 ) ' 0*1) 
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expected and observed intensity ratio before the first bakeout 
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Figure 40. Comparison between the expected (colored lines) and observed (black +) X-ray 
intensity ratio. The expected intensity ratio is calculated with the calibrated thicknesses of the 
filters and contaminant. 



where t o p - and t°'P- ±1 are the cumulative time when the filter located on and 
next to the optical path, respectively, and the subscripts of FPAF1 and FPAF2 
mean the FPAFs mounted on FW1 and FW2, respectively. On the basis of 
the estimated thicknesses of contaminant on the four thin filters (Al-mcsh, Al- 
poly, C-poly, and Ti-poly) in Phase 3 (see Table EJ) and the cumulative time 
t of filter usage shown in Figure I38l wc derive the accumulation rate of con- 
taminant on FPAFs as follows: x = 3220 A month -1 , y = 480 A month -1 , 
z = 2080 A month" 1 and w = 1000 A month" 1 . The results of x > y and z > w 
is likely, because the filter wheels are contained in a closed structure which has 
holes for the optical path and the contaminant may have come to the filters 
through such holes. 

So far, we derived coefficients of proportionality x, y, z, and w in Equa- 
tions ([30)1 and ([3Tj) using the final thicknesses of the contaminant for the four 
thin filters at the end of Phase 2 (beginning of Phase 3). If the thus-derived 
coefficients are correct, Equations ([3171) an d (EH should give correct contaminant 
thickness for each FPAF at any instant in Phase 2. Let us now check if this 
is the case. Figure [7] indicates the temporal evolution of contaminant thickness 
on the FPAFs from Equations ([30]) and ([3jj), while that for the CCD is shown 
in Figure [SJ Figure [30] shows comparison between observed intensity ratios in 
X-rays and those expected using the thus-derived contaminant thickness. Note 
the contaminant accumulated on both FPAFs and CCD in period A while only 
on the CCD in periods B-E. The close match between observed and expected 
intensity ratios strongly suggests that the estimate of contaminant thicknesses 
on the FPAFs with Equations ([3U|) and ([3Tj) are reasonable ones. We now have a 
reasonable estimate on the thickness of contaminants not only on the CCD but 
also on the FPAFs throughout the entire period since launch. 
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We note that, since not enough X-ray signal is acquired from the quiet Sun 
with the thicker filters, the contaminant thicknesses on the thicker filters have not 
been directly measured in contradiction to the thinner filters. Hence, the estimate 
of the thickness of contaminant on the thicker filters with Equations (f5U|) and 
(|3ip may be less accurate than the thinner filters in the sense that verification 
with X-ray data is not possible. Nevertheless, we do not see any reason for 
Equations ([50)1 and (|31[) not to be applicable to the thicker filters. Furthermore, 
the effect of contamination on the thicker filters is much smaller than on the 
thinner filters because only short-wavelength X-rays can transmit thicker filters 
and the effect of contaminant is small or negligible for such short wavelength 
X-rays. Hence, our results can be applied to quantitative analyses with XRT 
data. 

C. CCD Bakeout 

The XRT team decided that the 800 A of contaminant thickness should be 
the maximum acceptable contamination. The XRT team regularly performs 
CCD bakeout every three to four weeks to remove the contaminant from the 
CCD. Table [7] is the summary of CCD bakeout from the launch of Hinode (22 
September 2006) to the end of April 2009. 

D. Estimate of Data Number and Intensity Detected by XRT 

We explain the method to estimate the expected data numbers (DNs) when 
XRT observes the X-ray spectra from an X-ray generator and from the Sun, in 
order to compare them with the actually detected DNs at XRCF and on orbit, 
respectively. First of all, we define the meaning of DN and intensity in this article. 
The DN is the dimensionlcss value which is proportional to the total energy of 
the incident X-rays into XRT. Meanwhile, the intensity is the DN normalized 
by observation time and number of pixels [DN sec -1 pixel -1 ]. 

D.l. XRCF 

When incident X-rays into XRT have a photon-number spectrum of P (A) [cm~ 2 
sec" - 1 A" 1 ], the energy spectrum observed by XRT (£xrt (A) [eV sec 1 A 1 ]) 
is written as 



where h, c, and e are Planck's constant, speed of light, and elementary electric 
charge, respectively. This incident energy generates electron hole pairs on the 
CCD. One electron hole pair is generated by 3.65 eV. The data number DN 
detected by the XRT CCD is expressed as 



Exrt (A) = P (A) x A cB (A) x — x -, 



(32) 



DN 



t x 



/ 



( 



-Exrt (A) x 



3.65 



1 




(33) 
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Table 7. Bakeout history of XRT. 



No>) 


bake heater on 


heater off 


interval 


contam. rate 


spot 








[days]<W 


[A month- 1 ] (c) 




-3 


22 Sep 2006, 21:39 


16 Oct, 07:53 


0.1 




0.00 


— 2 


16 Oct 2006, 09:35 


17 Oct, 08:28 


0.1 




0.00 


_1 


17 Oct 2006, 10:13 


19 Oct, 08:12 


277.0 




0.00 


1 


23 Jul 2007, 09:09 


24 Jul, 08:10 


6.0 


- 


- 


2 


30 Jul 2007, 08:41 


3 Sep, 09:12 


148.0 


730 


2.61 


3 


29 Jan 2008, 08:42 


1 Feb, 10:01 


6.0 


1866 


5.48 


4 


7 Feb 2008, 08:50 


8 Feb, 07:52 


13.0 


716 


- 


5 


21 Feb 2008, 08:19 


22 Feb, 08:22 


13.0 


747 


5.25 


6 


6 Mar 2008, 08:18 


7 Mar, 02:20 


20.2 


613 


- 


7 


27 Mar 2008, 08:15 


28 Mar, 08:12 


20.0 


656 


- 


8 


17 Apr 2008, 09:10 


18 Apr, 10:00 


27.0 


569 


5.26 


9 


15 May 2008, 09:18 


16 May, 09:22 


13.0 


801 


5.25 


10 


29 May 2008, 09:50 


30 May, 10:00 


19.9 


506 


5.27 


11 


19 Jun 2008, 08:14 


20 Jun, 10:00 


20.0 


481 


5.28 


12 


10 Jul 2008, 09:47 


11 Jul, 04:00 


25.2 


432 


5.24 




Aug ZUUo, Uo. oo 


Aug, U4.UU 


ZZ.o 


OlU 


c: Q'3 
o.zo 


11 


28 Aug 2008, 10:05 


29 Aug, 04:00 


25.2 


514 




15 


23 Sep 2008, 09:53 


24 Sep, 04:00 


22.3 


522 


5.24 


16 


16 Oct 2008, 11:06 


16 Oct, 23:07 


20.4 


567 


5.24 


17 


6 Nov 2008, 09:29 


6 Nov, 21:29 


20.5 


507 


5.23 


18 


27 Nov 2008, 10:11 


27 Nov, 22:17 


20.4 


546 


5.23 


19 


18 Dec 2008, 08:33 


18 Dec, 20:39 


20.5 


544 


5.23 


20 


8 Jan 2009, 09:24 


8 Jan, 21:23 


20.5 


547 


5.23 


21 


29 Jan 2009, 08:29 


29 Jan, 21:23 


22.5 


487 


5.22 


22 


21 Feb 2009, 08:24 


22 Feb, 02:26 


18.3 


527 


5.22 


23 


12 Mar 2009, 09:19 


12 Mar, 21:20 


20.5 


499 


5.21 


24 


2 Apr 2009, 09:16 


2 Apr, 21:25 


20.5 


577 


5.23 


25 


23 Apr 2009, 09:13 


23 Apr, 21:14 






5.18 



( a ) Bakeouts from —3 to —1 were preformed before first light (23 October 2006). We count 
the bakeout number up with positive number after the first light. 

O Interval from this bakeout to next bakeout. 

' c ' Accumulation rate of contaminant on the CCD from this bakeout to the next bakeout. 
This rate is measured from G-band intensity oscillation. The details are described in 
Appendix IB. 2.2l 

(d) 'pjjg ra tj f contamination spot area to the full CCD area. The details of the spot are 
described in Appendix [Gl 
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where t is the exposure time in the unit of [sec], and G is the system gain of 
the CCD camera, which is 57.5 [e DN -1 ] in the XRT case QKano et al, 2008). 



Substituting Equation (|3"2"j) to (|3"3"|) . we obtain the equation to derive the DN 
from a photon- number spectrum P (A), 

DN = t,J (^(A)xA cff (A)x^X ex3 ] 5xG )dA. (34) 

D.2. On Orbit 

Using the effective area (see Figure^, we can estimate the intensity [/] detected 
by XRT, i.e. how much data number [DN] is detected by one pixel of XRT CCD 
in one second, when XRT observes the Sun whose photon-number spectrum is 
Pq (A) [cm -2 sec -1 sr -1 A -1 ]. The estimated / [DN sec -1 pixel -1 ] is given by 

^/(p^xSx^xf*-^)^, (35) 

where S [cm 2 ] is the solar area detected in one pixel of CCD. A c g and R are 
the effective area of XRT and distance between the Sun and XRT (« 1 AU), 
respectively. Then A e s/R 2 is the solid angle of XRT effective area in units of sr. 
The meaning of the other terms are the same as Equation (f3~4"|) . Considering the 
geometry of the telescope, the relation of 

w = f 2 (36) 

is given, where s and / are the area of CCD 1 pixel (13. 5 2 /im 2 , see lKano et al.\ 
POOS)) ) and focal length of XRT (2708 mm, see lGolub efHn([2W7j) ). respectively. 
Using this relation, Equation (|35[) can be rewritten as 

^ eff (A) he 1 

f 2 X A X e x 3.65 x G 



I = I ( P (A) x a x x - x _ oce w ^ ) dA. (37) 



Using Equation (|3T|) . we can derive the intensity observed by XRT. The CHI- 



ANTI atomic database (Dere et al., 19971 gives us a photon- number spectrum 
Aso© (A, T) [cm -2 sec -1 sr _1 A -1 ] emitted from an isothermal plasma at a given 
temperature of T and an unit column emission measure (CEM) of 1 cm -5 . The 
definition of CEM is 

CEM ee J n c x n H dl, (38) 

where n e , nn, and dl are the electron number density [cm -3 ], hydrogen number 
density [cm -3 ], and unit length along the line-of-sight [cm], respectively. We 
calculate P iso o (A, T) based on CHIANTI version 6.0.1 (jDere ' et alA IT9971 [20091 
which is the latest version, when we analyzed) with ionization equilibrium of chi- 
anti.ioneq ( Dere, 2007[) and abund ance of sun_coronal_ext. a bund ( |Feldman et al., 19921 1 
|Landi, Feldman, and Dere, 20020 . Replacing P Q (A) by P isoQ (X.T) in Equation (|3"7|T | 
we can calculate the so-called "temperature response" F (T) [DN sec -1 pixel -1 ] 
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of XRT, which is the intensity observed by XRT from a plasma at a given 
temperature of T and a unit column emission measure of 1 cm~ 5 as 

A eS (A) he 1 

f 2 X A X e x 3.65 x G 



F(T)= I ( P iso0 (A, T) x a x -2g-L x — x — ) dA. (39) 



Figure [10] shows F (T) for each X-ray analysis filter at the launch of Hinode 
and on-orbit. According to this figure, we can appreciate that the XRT has a 
superior ability to observe coronal plasmas over a wide temperature range from 
less than 1 MK to more than 10 MK. 

Using this temperature response F(T) of XRT, we can estimate the X-ray 
intensity / [DN sec -1 pixel -1 ] observed with XRT as follows: 

I = F(T) x CEM (40) 

for an isothermal corona at a temperature of T, and 

I = J F{T) x DEM (T) dT (41) 

for a multi-temperature corona with the differential emission measure DEM 
[cm~ 5 K _1 ], which is defined as 

DEM = - — ~j7^r ~ ^ ■ (42) 

The data number [DN] detected by XRT is given as 

DN = I xtxp, (43) 

where t and p are the observation time and the number of pixels used for ob- 
taining the DN, respectively. Using the volume emission measure [VEM: cm~ 3 ] 
defined as 

VEM = J n c x n H d^ dS = CEM x (S x p), (44) 

where dS is an unit area of observed solar region, DN is expressed as Equa- 
tion for an isothermal corona at a temperature of T. 



E. X-Ray Spectrum at XRCF 

We estimate spectra of X-ray beams at the XRCF. Figure [41] shows the con- 
figuration of the end-to-end test performed at the XRCF. The XRCF employs 
a target-impact type X-ray generator. Five X-ray lines that we used in this 
test are summarized in Table [8] with target sources, applied acceleration voltage 
for electrons, and blocking filters. The blocking filters, called E- and W-filters, 
are introduced to suppress longer and shorter wavelengths than the selected 
characteristic X-ray line, because lower levels of characteristic X-rays and con- 
tinuum bremsstrahlung X-rays are also emitted. The X-rays illuminate XRT in 
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X-ray 




Figure 41. Configuration of the end-to-end test at XRCF. 



Table 8. Characteristic X-ray lines used in the end-to-end test at XRCF. 



line 


energy 


wavelength 


source 


voltage 


E-filter 


W-filter 


Mo-L 


2.29 keV 


5.41 A 


Mo mono 


7.0 kV 


Lexan( a ' 


Mo 3.0^m 


Al-K 


1.49 keV 


8.34 A 


Al mono 


7.0 kV 


Lexan 


Al 18.7^m 


Cu-L 


0.930 keV 


13.3 A 


Cu mono 


1.45 kV 


Cu 2.8/^m 


Lexan 


O-K 


0.525 keV 


23.6 A 


O on Al 


4.0 kV 


Cr 0.7fim 


Lexan 


C-K 


0.277 keV 


44.7 A 


C mono 


5.0 kV 


Lexan 


Parylene 10.0/^m 



W This Lexan is the product of Luxel R/N 6370. 



an almost parallel beam, because the small-sized X-ray source is located far from 
XRT, namely ~ 500 m. A flow proportional counter (FPC) was located beside 
the XRT to monitor intensity and spectral shape of the X-ray beam. 

Here, using the FPC data, we derive the incident X-ray spectrum into XRT. 
The FPC detected two humps in each characteristic X-ray line data (see Fig- 
ure 02) ■ Even with the blocking filters, the continuum bremsstrahlung emis- 
sion component is not completely removed. The main hump comes from the 
characteristic X-ray line and the other from the continuum bremsstrahlung. 

To understand the incident X-ray spectrum into the XRT, we need to estimate 
the spectrum of continuum bremsstrahlung also. On the basis of the thick target 
emission model, the continuum bremsstrahlung profile can be described as 

/thick(e) = 4dM j F ( E o)J Ea B (e,E)dEdE , (45) 



a B (e,E) = 1.58 x RT 24 -hi 

£ 

QSakao, 1994[ ) , where e is the photon energy, E the electron energy, E the initial 
electron energy, F the energy distribution function of electrons injected to the 
target [electrons sec -1 keV -1 ], ctb the Bethe-Heitler cross section ( | Jackson, 1962[ ),| 
K = 27re 4 lnA, In A the Coulomb logarithm ( |Spitzer, 1962[ ), and R the distance 
between the X-ray source and the FPC. In our case, because single-energy 
electrons at an energy of Eqq were injected, the distribution function can be 
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Figure 42. X-ray spectra detected by the flow proportional counter (FPC). The black marks 
(asterisks) are the data points detected by FPC. The orange line shows the estimated spectrum, 
which consists of characteristic X-ray line (red line) and continuum bremsstrahlung (blue line). 



expressed as: 

F(E )=F ao 6(E -E ao ), (47) 

where Fqq is the number of injected electrons per second. The dashed line in 
Figure [43] is the estimated continuum bremsstrahlung spectrum emitted at the 
acceleration voltage of 7 kV with Equations (|4"5]l - (|4"T|) . Using the transmission of 
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continuum bremsstrahlung in Mo-L case 

-i 1 1 1 1 1 1 1 1 1 1 1 1 1 r - 




energy [keV] 

Figure 43. Continuum bremsstrahlung spectrum. The dashed line is the continuum spectrum 
emitted by 7 keV electrons. The dotted line indicates the incident continuum spectrum into 
FPC. The solid line shows the continuum spectrum detected by FPC with its energy resolution. 

Table 9. Best-fit photon number fluxes of characteristic X-ray 
lines and continua. 



line 


energy 


voltage 


characteristic line 


continuum 




[keV] 


[kV] 


[cm -2 sec -1 ] 


[cm -2 sec -1 ] 


Mo-L 


2.29 


7.0 


80.7 (63.6%) 0) 


46.1 (35.4%) 


Al-K 


1.49 


7.0 


77.4 (90.8%) 


7.8 (9.2%) 


Cu-L 


0.930 


1.45 


50.8 (36.1%) 


89.9 (63.9%) 


O-K 


0.525 


4.0 


2.2 (68.7%) 


1.0 (31.3%) 


C-K 


0.277 


5.0 


7.8 (30.9%) 


17.4 (69.1%) 



( a ) The percentage in parenthesis shows the ratio of each component to the total photon flux. 



the E- and W-filtcrs and the efficiency of FPC, the incident continuum spectrum 
into FPC is derived as shown by the dotted line in Figure S3] Considering the 
energy resolution of FPC, which can be written as 



(AE) 



FWHM 



= aE 1 ' 2 



(48) 



QCharles and Cooke, 1968| , where a is a constant of proportionality, the contin- 
uum spectrum detected by the FPC is simulated as shown by a solid curve in 
Figure |H 

The actual spectrum detected by the FPC is fit by this simulated FPC 
continuum spectrum and the characteristic X-ray line (Figure |42]) . The fitting- 
parameters are the constant a in Equation (|4"8"1) . the number of injected electrons 
per second Fqq , and the strength of the characteristic X-ray line represented by 
a Gaussian shape. 
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Figure 44. Stainless mesh of Al-mesh filter taken by the filter manufacturer with a high 
precision digital microscope. This mesh supports the thin Al film of Al-mcsh filter. 



iijijSjjjiji b 




Figure 45. Concept images of variation in "annulus transmission". The gray area indicates 
the mesh of Al-mesh filter. The black annulus shows the shape of X-rays imaged on the filter. 
The left and right panels shows the difference of annulus transmission caused by the different 
relative position between mesh pattern and annular X-rays. The aspect ratio between the 
black annulus and gray mesh in this figure is different from the actual aspect ratio between 
the annular X-rays image and stainless mesh on Al-mesh filter. 



Using the best-fit results in Figure WZ[ we obtained the incident spectrum into 
XRT (sec Figure ITO)) . The errors shown by gray areas in Figure \W\ arc estimated 
from the photon noise and the energy range of FPC data points. For example, 
the error bars in 0-K and C-K line are larger than for the other lines, because the 
total photon number is small in the O-K line case, and because the energy range 
of the C-K line does not completely cover the continuum. The photon-number 
fluxes of characteristic X-ray lines and continua are summarized in Table [HI The 
contribution of continuum is not negligible even for the Al-K line. 



F. Variation in Annulus Transmission 

The support of the Al-mesh filter is a stainless steel mesh. Figure PHI is a picture 
of the mesh taken by the filter manufacturer with a high-precision digital micro- 
scope. According to the manufacturer, the width of the mesh wire and the open 
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Figure 46. Variation in "annulus transmission". The left panel is the mesh cell. The black 
and white regions indicate the mesh wire and open area, respectively. The right panel shows the 
annulus transmission where the center of annular X-rays imaged is located at the corresponding 
left panel. The white and black mean high and low transmission, respectively. 
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Figure 47. Histogram of annulus transmission based on the right panel of Figure I46I The 
dotted line indicates the average of annulus transmission: 77.1%. 



area are 38 /im and 330 /im, respectively. On the basis of this measurement, the 
open area occupies 77% of whole area. We adopt this measured value of 77% for 
the geometrical open area. 

Here we define "annulus transmission" as the fraction of the X-rays passed 
though the open area of the mesh to the incident X-rays into the mesh. Note 
that the X-rays pass through the open area without any loss of intensity, while 
at the wires, the X-rays are completely blocked. Hence, the average of annulus 
transmission should be consistent with the geometrical open area of the mesh. 

The annulus transmission is not uniform and depends on the relative position 
of the shape of X-rays on the filter with respect to the mesh pattern. However, 
we do not know this relative position for the flight XRT. Instead, we derive 
distribution of the annulus transmission for all possible relative positions, take 
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the average of the distribution as the annulus transmission of the Al-mesh filter, 
and employ the deviation of the distribution as that for the annulus transmission. 
The procedure is detailed in the following. 

The shape of X-rays imaged on the FPAFs is an annulus, because the shape 
of the XRT aperture is an annulus (see Figure [3]). The radius and width of 
annular X-ray image are about 4.1 mm and 7.6 fim on the FPAFs, respectively, 
because the radius and width of XRT annular aperture, focal length, and distance 
between focus position and FPAF are about 170 mm, 0.32 mm, 2700 mm, and 
65 mm, respectively. Figure 05] shows some schematic examples of the relation- 
ship between the annular X-ray image and the mesh. The gray area indicates 
the mesh of Al-mesh filter and the black ring shows the X-ray annulus on the 
filter. The variation in annulus transmission is shown in Figure 1461 The left 
panel shows a cell (unit opening) of the mesh. The width and height of this cell 
are both 368 fj,xn on the basis of manufacturer's measurement. The black and 
white regions indicate mesh wire and open area, respectively. The right panel 
shows the annulus transmission in gray scale at the position where the center 
of the annular X-ray image is located at the corresponding left panel. Figure l47l 
shows distribution of the annulus transmission over the varying center position 
of the X-ray annulus with respect to the mesh pattern (This is the histogram of 
calculated annulus transmission in the right panel of Figure I46[) . The average of 
annulus transmission is consistent with the geometrical open area of the mesh. 
However, the annulus transmission has a scatter with the standard deviation a 
of 2.3%. This a is also considered in the calibration of Al-mesh filter thickness. 

The uncertainty in the annulus transmission originated from the known rela- 
tive position of the annulus center with respect to the mesh pattern can thus be 
evaluated in this way. In this evaluation, we assumed that the shape of the X-ray 
image on the filter is a perfect annulus. But, in reality, the shape corresponds to 
the colored area in Figure [3] and is not symmetrical. Hence, we should consider 
not only the relative position but also the relative angle between the patterns 
along the image annulus and the direction of the mesh wires, which is also not 
known. In order to assess the amount of additional uncertainty caused by this, 
we calculated the scatter of annulus transmission due to the uncertainty of the 
relative position and angle. The standard deviation a' of this scatter (imperfect 
annulus case) turned out to be almost the same as a = 2.3% (perfect annulus 
case). Hence, we claim our treating the image shape as a perfect annulus is 
sufficient for deriving the average fraction and a of the annulus transmission for 
the X-ray image on the Al-mesh filter. 

G. Spot-Shaped Contaminants on the CCD 

After the first bakeout, spot-like patterns appeared in both X-ray and G-band 
data as shown in Figure l48l Subsequent analyses indicated that the estimated 
thickness of contaminant on the CCD was more than 3000 A before the first 
bakeout (see Figure [5]). The patterns most likely originate from the condensa- 
tion of contaminant and are triggered by a bakeout with the thickness of the 
contaminants accumulated on the CCD exceeding a certain threshold value (the 
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spots in X-ray image 



spots in G-band image 







Figure 48. Spot-shaped contaminants observed in X-rays and G-band. The spots in the 
X-ray data are dark, because the spots absorb the X-rays. Meanwhile, the spots in the G-band 
data are bright, since the spots work as anti-reflection coatings. 



value itself is not known). These spots were not removed at all even with the 
second bakeout that was performed after the first bakeout and lasted for a month. 
The ratio of spot area to full CCD area was 2.6% at this point (see Table [7]). 
As it took five months before we agreed, as the instrument team, to proceed 
with performing the next (third) bakeout, the resultant thickness of the con- 
taminants again exceeded 3000 A, reaching even beyond 3500 A (sec Figure H2). 
This resulted in the creation of additional contamination spots across the CCD, 
with an increased spot area ratio of 5.2%. Note, however, this additional increase 
was rather an intentional one; We took the option of keeping a low-temperature 
diagnostic capability by removing the accumulated contaminants while accepting 
an increase in the spot area. The resultant spot distribution is indicated by the 
black areas in FigureHHl After the third CCD bakeout, we regularly perform the 
bakeout every three to four weeks. As we guessed, no more spots were created 
after the third bakeout as shown in Table [7] The ratio of spot area to full-CCD 
area has kept to be 5.2%. 

We investigate the relation between the thickness of some spots measured 
with their Newton-ring patterns in G-band and the absorbed X-ray intensities 
by such spots, and conclude that the material of spot contaminant is most likely 
the same material as the laminar contaminant. This result will be presented as 
a separate paper with more detailed analysis. 



H. Transmission of Rays Through Three Layers 

In Appendix IB. 2. 11 in order to measure the contaminant thickness on the CCD, 
we used the transmission of visible light into the CCD through contamination 
layers, Equation (|22[) . Here, we derive this equation. 
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spot map after the third bakeout 




Figure 49. Spot map after the third bakeout. The ratio of spot area to full CCD area is 
about 5%. 



The reflected and transmitted amplitudes [Ar and At] of rays at a boundary 
of two layers whose refractive indexes are rii and n 2 are written as 

Ah = A f^l) , (49) 
\ni + n 2 J 

A T = A (-^-) , (50) 

where A is the amplitude of incident rays from the n\ layer to the boundary. 

We consider transmission of three layers as shown in Figure 1251 Here, we 
define A as the amplitude of incident rays in vacuum (m layer), and Ak the 
complex amplitude of rays which are reflected 2k times at the borders of layers 
(rii — 712 border and n 2 -ri3 border) and transmitted into the CCD (713 layer). 
Note the amplitudes of rays shown by white and black arrows in Figure [23 are 
Aq and Ai. Ak is expressed as 

Ak = Ae MA, ( ( r ( ( y (51) 

V n-i +n 2 J \n 2 +n 3 / I V n 2 + n z I \n 2 + m J I 
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where Aip (= &2 x 2d) gives the phase difference between Ak and Ak+i- Hence, 
the total amplitude [At] of incident rays into the CCD is derived as 



2ni \ / 2n 



k=0 



ri\ + ri2 ) \n>2 + 113 



1 _ e iA V { ( n 2~ n A f n 2 - ni 
\n 2 +n 3 y Vti 2 +ni 



(52) 



Using the relation between amplitude A and intensity I in a layer of refractive 
index n given as 

X = nAA*, (53) 

where * denotes complex conjugate, the intensity / in vacuum (m layer in 
Figure 123)) is 



I = mAA*, (54) 
and the total intensity I t detected by the CCD is 

I t =n 3 A t Al. (55) 

Equation (|22|) is derived from Equations (|52|) . (|54l) . and ([55]l, 



I. Oxidization of Metal 

When the molecular formula of pure metal and oxidized metal are X and XaOb, 
and the molecular weight of X and O are Mx and Mo(= 16.00), respectively, 
the reduced thickness [Ac? pure ] of pure metal by the creation of oxidized metal 
with a thickness of <i x is derived to be 

Adp Uro = d ox x — -, 56) 

Mx x yl + Mo x £? 

based on the conservation of the number of metal atoms. The information about 
materials used in our calibration is summarized in Table [4] 



J. Meaning of Photon Noise Calculated with Filter-Ratio 
Temperature 

The photon noise can be derived with Equation (|13p . For the Condition (B) in 
Section I5~2l we investigated how much photon noise is expected in the observed 
data with the single temperature of the filter-ratio temperature. However, the 
actual solar corona may have a multi-temperature structures. Hence, here we 
discuss the meaning of photon noise estimate with the filter-ratio temperature. 
Figure [T7] shows the correlation between the fractional photon noise derived with 
DEM [a£§ M /DN] and with the filter-ratio temperature [a^/DN] for four 
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Correlation between the fractional photon noise derived with DEM [c^^ M /DN] 



and with the filter-ratio temperature [cr^^° /DN] for four DEMs shown in Table[3] This plot 



is the case when 1 month passed after the CCD bakcout, i.e. 800 
CCD. The dashed line indicates the positions whe 
number of analyzed filters is shown in parenthesis. 



of contaminant on the 

CCD. The dashed line indicates the positions where a™^°/DN is equal to a°^ M /DN. The 



DEMs (regions) summarized in Table [3] erff^ 10 is the photon noise calculated 
with the filter-ratio temperature [T rat i ] and Equation (I13|) . crgg M is the photon 
noise derived with the DEM as following steps: i) We calculate the photon noise 
[cdn] for each temperature bin of DEM with Equation (fT"3"]) . ii) We derive the net 
photon noise of step i) as cr^ M = \/J2 °dn- ^ ne ct dn° an d °dn M ^ s calculated 
for all filters which are included in the suitable filter pairs in Figure [H] with the 
typical values of coronal structures in Table [31 On the basis of this figure, we 
can say that crjy^° well matches crp§ M which is the actual photon noise. 

This equality can be understood as following: Since the amount of plasma 
at the DEM peak is dominant, its photon noise mainly affects ct§n M - ^ n the 
temperature range where the DEM peak is located, the conversion factor [KW] 
from observed DN to ^/<jdn (see Equation [T3")) is almost constant for the filters 
which are included in the suitable filter pairs (see Figures [T3l and [14]) . For exam- 
ple, in a temperature range of 3-4 MK where the active region has the DEM 
peak (see Figure [16]), for Al-poly, Ti-poly, and med-Be filters varies within 
factors of 1.07, 1.11, and 1.02, respectively (see the bottom panel in Figure IT2"]) . 
Hence, cr^ M well matches a™$°, since dj^ io is derived from at T rat i , and 
Iratio is close to the DEM peak temperature (see Section 153)) . We note that, 
for temperature ranges other than the DEM peak, varies with a factor of 
typically 2-4. However, since i) the photon noise is a function of VA® (see 
Equation I13[) and furthermore, ii) in this temperature range, the contribution 
of photon noise to is smaller than at the DEM peak, the variation of 

in these temperature ranges does not significantly affect <r|]™. 
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K. Notes 

K.l. Contaminant Material 

In Appendix IB. 3. 11 wc identified the contaminant material as a long-chain or- 
ganic compound whose chemical composition, density, and refractive index are 
similar to those of DEHP, although DEHP is not used in the XRT. This means 
that the actual contaminant material is narrowed down to the material whose 
refractive index and density are close to 1.5 and 1 g cm~ 3 , respectively. 

Next let us discuss the chemical composition of the actual contaminant mate- 
rial. Since the atomic number of silicon is about two times larger than carbon and 
oxygen, i.e. the cross-section of silicon to X-rays is certainly larger, the existence 
of silicon affects on the X-ray transmission certainly. In fact, when we assume 
that the contaminant contains some silicon atoms, the discrepancy between the 
observed and expected X-ray intensity ratios rises as shown in Figure l33l Hence, 
we conclude that there are no silicon atoms in the contaminant material. 

On the other hand, since the atomic numbers of carbon and oxygen (also 
nitrogen) are comparable and that of hydrogen is much smaller than carbon and 
oxygen, i.e. hydrogen is much more transparent to X-rays, some slight differences 
in numbers of such atoms among a long-chain organic compound do not make 
any significant difference in the transmission of X-rays. Hence, although we 
cannot exactly identify the chemical composition of the contaminant material, 
it is acceptable for the calibration of the XRT that we tentatively employed 
DEHP as the contaminant material. We note that the identified property of 
the contaminant, where chemical composition is a long-chain organic compound 
without silicon, refractive index is 1.5, and density is sa 1 g cm -3 , is a common 
property of materials which are widely used for satellites, and which are well 
known as possible source of contamination. 

K.2. G-Band Method for CCD Contamination Analysis 

In Appendix IB . 2.11 wc explained the method to measure the thickness of contam- 
inant on the CCD using G-band data. The G-band intensity oscillation shown 
in Figure HH is essentially expressed by Equation (j2"5)l . However, there are two 
differences between the calculated G-band intensity enhancement with the anti- 
rcflcction effect of Equation (|23|) and the observed one that is characterized by 
the empirical model given by Equation (|26|) . 

The first one is the refractive indexes of the contaminant. In this article, we 
obtained the refractive index of contaminant with following two methods: 

i) In Appendix lB.2.2l from the observed amplitude of intensity enhancement 
in G-band, the refractive index of contaminant is estimated as « 1.4. 

it) In Appendix IB. 3. 11 the decrease in X-ray intensities caused by the accu- 
mulation of contaminant can be well explained by adopting a well-known 
material whose refractive index and density of contaminant are 1.5 and 
0.986 g cm -3 , respectively, as the contaminant material. 
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This difference may be explained with the possible hypothesis where the density 
of actually-accumulated contaminant is lower than the contaminant in its usual 
state (as, e.g., products), because the contaminant material is an organic com- 
pound and its accumulation was formed under vacuum deposition process but 
with much slower rate than standard ones in laboratories, and the accumulated 
contaminant most likely consists of sparsely-structured molecules (rather than 
the dense crystalline structure). Hereafter, we call this state as sparse state. The 
same trend as this hypothesis is seen in the case of C-poly filter as described in 
the last paragraph of Appendix IA.2I Generally, the refractive index [n S p arse ] and 
density [p spa rso] of material in sparse state are smaller than those (n usua i and 
Pusuai) m the usual state. If the above hypothesis holds, the method i) gives the 
refractive index of contaminant in sparse state. On the other hand, in method ii), 
we identified the refractive index and density in usual state from among the 
possible candidates of contaminant (see Table [5]). Hence, the difference in the 
refractive indexes in methods i) and ii) is caused by the different state of the 
contaminant. 

Next, we consider the effect of this difference on the calibration results. On the 
basis of Equation (|25p . the smaller refractive index gives the thicker contaminant. 
Meanwhile, the smaller density gives the longer attenuation length [i att ] in Equa- 
tion (|20[) . Hence, the X-ray transmission calculated with n sparse and p sparS e would 
be close to the transmission with n usua i and p U suai- We note that even if the actual 
state of contaminant is not sparse, i.e. even if the density of actual contaminant 
is the same as the usual state, the difference in refractive indexes (1.4 and 1.5) 
is only 7%, and then the difference in the estimated thickness of contaminant 
on the CCD is also 7%. Since the maximum thickness of contaminant is about 
3600 A (see Figure [S]), the maximum error in the estimate is about 250 A. This 
error is smaller than the accumulated thickness of contaminant on the thinner 
FPAFs (see Figure [7]) , and is negligibly small against the metal thicknesses of 
the thicker filters. Hence, we conclude that, whether the hypothesis holds or 
not, the difference in refractive indexes does not significantly affect the results 
of XRT calibration, and adopt n usua i = 1.5 and /0 U suai = 0.986 g cm -3 for the 
calibration. 

The second one is the decay of oscillation amplitude. We consider that this 
is caused by the inhomogeneous accumulation rate of the contaminant. In Ap- 
pendix IB. 2. 31 we demonstrated the accumulation rate of contaminant at the 
CCD center is larger than at the edges of the CCD. This effect might appear 
in the 512 x 512 pixels data which is used to make the G-band intensity plot 
(Figure |2"4")) . The difference of accumulation rate will result in an asynchronous 
phase in the intensity oscillation across the 512 x 512 pixels. The G-band inten- 
sity plot is the average of such different oscillations which are sine curves with 
different phases. Because the phase difference becomes larger and larger as time 
passes, the oscillation amplitude of the G-band intensity enhancement decays. 

On the basis of Equation (|2"5j) . the thickness of contaminant on the CCD is 
derived from only the refractive index of the contaminant and the period of the 
intensity oscillation. Hence, we expect that the refractive index of ss 1.5 and the 
oscillation period in Figure [53] give the average thickness of contaminant on the 
CCD area of 512 x 512 pixels, which is reliable enough for the XRT calibration. 
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K.3. Accumulation Profile of Contaminant on the CCD 



When we estimated the thickness of the contaminant accumulated on the CCD 
with the visible light intensity profile (Appendix IB . 2 . 2"|) . we assumed a constant 
accumulation rate for the contaminant between two successive bakeouts. How- 
ever, actually, we see rapid accumulation of contaminant right after the CCD 
bakcout while for the rest of each period the rate is almost constant (see right 
panel of Figure |2"9"| . Nevertheless, the assumption of constant accumulation for 
the entire period between two successive bakeouts is a valid one for the calcu- 
lation of the XRT effective area, because the contaminants accumulated on the 
FPAFs are much thicker than the difference between the actual and estimated 
thickness with the above assumption, i.e. such difference is negligible. 

K.4. DEM Model Used for the Analysis 

In Appendix IB. 31 we calibrated the contaminant thickness on the FPAFs with 
the observed X-ray data and DEM in the quiet Sun. When we performed the 
calibrations described in this article, the DEM model in the quiet Sun derived 
bv lBrooks and Warrenl (|2006|) was the latest result available. Hence, we adopted 
their DEM model for our analysis. IBrooks and Warrenl (|2006p derived the DEM 
of quiet Sun from a data set taken with SOHO/E IT, SOHO /CDS, and TRACE 
on 1 May 1998, while more recently fBrooks et aZ.I(|2009p analyzed the DEM with 
45 data sets observed with Hinode /EIS in the period from January to April 2007. 
The profiles of DEMs studied in these two papers are very similar up to at least 
logT rj 6.2 K in all cases. We confirmed that the calibrated contaminant 
thicknesses on the FPAFs with the DEMs from Brooks and Warren ( 2006) and 
IBrooks et dl\ (P0"0"9"|) are consistent with each other. 



K.5. Source Location of the Contaminant 



We found that there is remarkable similarity between the observed spatial distri- 
bution of the contaminant across the CCD and the simulated result by Urayama et aZ.[ 



(2008) for the case where the contaminant was assumed to come from the 
direction in front of the CCD (see Figure [30]). Meanwhile, we found that the 
contamination is triggered by the operational heater which warms the rear end of 
telescope tube up to about 20°C. On the basis of these circumstantial evidences, 
we suspect that the contaminant most likely originates from somewhere in front 
of the CCD, in the telescope tube. 

K.6. Future plan for the calibration 

As we mentioned in Appendix IA.2I we consider that the calibration of thick 
filters (med-Be, thick- Al and thick-Be filters) with ground-based test data is not 
ideal, though the difference between calibrated and actual X-ray transmission 
of them should be small. For further calibration of thick filters, we need to take 
enough data sets where XRT observes intense X-ray sources, e.g., active regions 
and flares, with thick filters. Calibration with on-orbit observation data is our 
future work. 
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Also, in order to supplement our calibration, we plan to perform the following 
cross-calibrations: i) between Hinode/XRT and Hinode/ElS, and ii) between 
XRT and GOES 13 /SXI, which is the grazing incidence X-ray telescope like 
XRT. 

Acknowledgements The authors thank members of the XRT team for useful discussions 
and comments. We acknowledge D. Brooks and H. Warren for variable information on quiet Sun 
and active region DEMs. Hinode is a Japanese mission developed and launched by ISAS/JAXA, 
collaborating with NAOJ as a domestic partner, NASA and STFC (UK) as international 
partners. Scientific operation of the Hinode mission is conducted by the Hinode science team 
organized at ISAS/JAXA. This team mainly consists of scientists from institutes in the partner 
countries. Support for the post-launch operation is provided by JAXA and NAOJ (Japan), 
STFC (UK), NASA, ESA, and NSC (Norway). CHIANTI is a collaborative project involving 
the NRL (USA), the Universities of Florence (Italy) and Cambridge (UK), and George Mason 
University (USA). We wish to express our sincere gratitude to the late Takeo Kosugi, former 
project manager of Hinode (SOLAR-B) at ISAS, who passed away in November 2006. Without 
his leadership in the development of Hinode, this mission would have never been realized. We 
express our sincere gratitude to those who participated in, or supported, XRT end-to-end 
calibration measurement at MSFC XRCF that was carried out for two weeks in May-June 
2005. First of all, we are very much indebted to an extremely-talented team of people at 
XRCF led by C. Reily and J. McCracken, including R. Siler, E. Wright, J. Carpenter, J. 
Keegley and G. Zirnstein, and also M. Baker, H. Haight, B. Hale, T. Hill, B. Hoghe, D. Javins, 
J. Norwood, H. Rutlcdgc, G. St. John, J. Tucker, and D. Watson. Without their continuous 
and enormous support including preparation of LN2-cooled cold plate facing XRT radiators 
in the cryogenic chamber of XRCF, from early preparation phase of the XRCF experiment 
and throughout the experiment period, the filter calibration at XRCF would have been totally 
impossible. We are also grateful to people from SAO who supported the XRCF measurement; 
G. Austin, W. Podgorski, E. Dennis, J. Chappel, D. Caldwell, W. Martell, M. Harris, M. 
Cosmo, D. Weaver, S. Park and T. Kent who participated in the experiment, and P. Cheimcts, 
J. Bookbinder, J. Boczcnowski. Also A. Sabbag of Naval Research Laboratory is appreciated 
for his support during prc-shipment instrument check performed at SAO. K. Kumagai and M. 
Tamura of NAOJ, K. Yaji of Rikkyo U., and K. Kobayashi of MSFC are greatly appreciated 
for preparing for, and participating in, the XRCF experiment. Also, people from MSFC, L. 
Hill, J. Owens, B. Cobb, D. Coleman, R. Jayroe, T. Perrin, D. Schultz, A. Sterling and C. 
Talley are sincerely acknowledged for their various support throughout our stay in Huntsvillc 
and also for our transportation to and from XRCF. Finally, the authors grateful to the referees 
and editor for their detailed and careful reviews. 

References 

Acton, L.W., Weston, D.C., Bruner, M.E.: 1999, J. Geophys. Res., 104, 14827. 
Brooks, D.H., Warren, H.P.: 2006, Astrophys. J. Supp. Ser., 164, 202 

Brooks, D.H., Warren, H.P., Williams, D.R., Watanabe, T.: 2009, Astrophys. J., 705, 1522 
Charles, M.W., Cooke, B.A.: 1968, Nucl. Instr. and Meth., 61, 31. 

Collura, A., Barbera, M., Inzerillo, G, Mirabello, F., Sciortino, S., Serio, S.: 1994, SPIE, 2280, 
206. 

Dere, K.P., Cook, J.W.: 1979, Astrophys. J., 229, 772. 



SOLA: Narukage_etal_R3.tex; 15 November 2010; 1:16; p. 75 



N. Narukage et al. 



Dere, K.P., Landi, E., Mason, H.E., Monsignori Fossi, B.C., Young, P.R.: 1997, Astron. 

Astrophys. Suppl., 125, 149. 
Dere, K.P.: 2007, Astron. Astrophys., 466, 771. 

Dere, K.P., Landi, E., Young, P.R., Del Zanna, G., Landini, M., Mason, H.E.: 2009, Astron. 
Astrophys., 498, 915. 

Fcldman, U., Mandelbaum, P., Seely, J.F., Doschek, G.A., Gursky, H.: 1992, Astrophys. J. 

Supp. Ser., 81, 387. 
Freeland, S.L., Handy, B.N.: 1998, Solar Phys., 182, 497. 

Golub, L., Deluca, E., Austin, G., Bookbinder, J., Caldwell, D., Cheimets, P., et al: 2007, 
Solar Phys., 243, 63. 

Hara, H., Tsuneta, S., Lemen, J.R., Acton, L.W., McTiernan, J.M.: 1992, Pub. Astron. Soc. 
Japan, 44, L135. 

Henke B.L., Gullikson E.M., Davis J.C.: 1993, Atomic Data Nuclear Data Tab., 54, 181. 
Jackson, J.D.: 1962, Classical Electrodynamics, J.Wiley and Sons, New York. 
Kano, R., Tsuneta, S.: 1995, Astrophys. J., 454, 934. 

Kano, R., Sakao, T., Hara, H., Tsuneta, S., Matsuzaki, K., Kumagai, K., et al: 2008, Solar 
Phys., 249, 263. 

Kosugi, T., Matsuzaki, K., Sakao, T., Shimizu, T., Sone, Y., Tachikawa, S., et al: 2007, Solar 
Phys., 243, 3. 

Landi, E., Feldman, U., Dere, K.P.: 2002, Astrophys. J. Letters, 139, 281. 
Martens, P.C.H., Cirtain, J.W., Schmelz, J.T.: 2002, Astrophys. J. Letters, 577, 115. 
Mcwc, R., Gronenschild, E. H. B. M., van den Oord, G. H. J.: 1985, Astron. Astrophys. Suppl, 
62, 197. 

Mewe, R., Lemen, J. R., van den Oord, G. H. J.: 1985, Astron. Astrophys. Suppl, 65, 511. 
Ogawara, Y., Takano, T., Kato, T., Kosugi, T., Tsuneta, S., Watanabe, T., et al: 1991, Solar 

Phys., 136, 1. 
Osantowski, J.F.: 1982, Proc. SPIE, 338, 80 

Powell, F.R., Vedder, P.W., Lindblom, J.F., Powell, S.F.: 1990, Optical Engineering, 29, 614. 
Sakao, T.: 1994, Ph.D. Thesis., The University of Tokyo. 

Schmelz, J.T., Kashyap, V.L., Saar, S.H., Dennis, B.R., Grigis, P.C., Lin, L., et al: 2009a, 

Astrophys. J., 704, 863. 
Schmelz, J.T., Saar, S.H., DeLuca, E.E., Golub, L., Kashyap, V.L., Weber, M.A., et al: 2009b, 

Astrophys. J. Letters, 693, 131. 
Schmelz, J.T., Saar, S.H., Nasraoui, K., Kashyap, V.L., Weber, M.A., DeLuca, E.E., et al: 

2010, Astrophys. J., 723, 1180. 
Shimizu, T., Katsukawa, Y., Matsuzaki, K., Ichimoto, K., Kano, R., Deluca, E.E., et al: 2007, 

Pub. Astron. Soc. Japan, 59, 845. 
Spitzer, L. Jr.: 1962, Physical of Fully Ionized Gases, 2nd ed., Interscience, New York. 
Tsuneta, S., Acton, L., Bruner, M., Lemen, J., Brown, W., Caravalho, R., et al: 1991, Solar 

Phys., 136, 37. 

Urayama, F., Bando, T., Kano, R., Hara, H., Narukage, N., Sakao, T.: 2008, Japan Soc. 

Aeronaut. Space Science, 56, 536. 
Vaiana, G.S., Krieger, A.S., Timothy, A.F.: 1973, Solar Phys., 32, 81. 
Vernazza, J.E., Reeves, E.M.: 1978, Astrophys. J. Supp. Ser., 37, 485. 
Warren, H.P., Brooks, D.H., Winebarger, A.R.: 2010, arXii;:1009.5976vl. 



SOLA: Narukage_etal_R3.tex; 15 November 2010; 1:16; p. 76 



